site 
iit 
Dall 


Gx BRIS 
UNIOMASITATIS 
ANBERTAEASTS 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Razdan1982 


DERY, Edmonton 


THE UNIVERSITY OF ALBERTA 
RELEASE FORM 


NAME OF AUTHOR RAJENDER RAZDAN 
TITLE OF THESIS AN INVESTIGATION OF MULTILEVEL PAM AND PWM FOR 
FIBER-OPTIC COMMUNICATION SYSTEMS 
DEGREE FOR WHICH THESIS WAS PRESENTED MASTER OF SCIENCE 
YEAR THIS DEGREE GRANTED SPRING 1982 
Permission is hereby granted to THE UNIVERSITY OF ALBERTA LIBRARY to 
reproduce single copies of this thesis and to lend or sell such copies for private, 
scholarly or scientific research purposes only. 
The author reserves other publication rights, and neither the thesis nor 
extensive extracts from it may be printed or otherwise reproduced without the 


author's written permission. 


THE UNIVERSITY OF ALBERTA 


AN INVESTIGATION OF MULTILEVEL PAM AND PWM FOR FIBER-OPTIC 
COMMUNICATION SYSTEMS 
by 
RAJENDER RAZDAN 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 
ORMASTER OF ISCIENEE 


ELECTRICAL ENGINEERING 


EDMONTON, ALBERTA 
SPRING 1982 


=» 


- : 7 
oj ey Aye NA RAT te) (GrTARNT aga wid 
VST OTA LAUDS 
— 


Sela gaoees, 47 


’ 
: a s @ 
e aie | ral 
uf 
5 MwA hae, ta! Sr aah a yteloa ge er catia 


SPlchettd> ph SRE Sut 4g Tae. Jal 


(aid 
=SS RO PTS al AC 


eos MST 


AT BA sgn 


sed) Quntiae . 
ae | 


THE UNIVERSITY OF ALBERTA 
FACULTY OF GRADUATE STUDIES AND RESEARCH 


The undersigned certify that they have read, and recommend to the Faculty 
of Graduate Studies and Research, for acceptance, a thesis entitled AN INVESTIGATION 
OF MULTILEVEL PAM AND PWM FOR FIBER-OPTIC COMMUNICATION SYSTEMS 
submitted by RAJENDER RAZDAN in partial fulfilment of the requirements for the degree 
of MASTER OF SCIENCE. 


ABSTRACT 

Multilevel transmission is advantageous when compared with binary transmission, 
since, for a channel of fixed bandwidth, more information can be transmitted. In this 
thesis, a theoretical investigation of the noise present in multilevel PAM and PWM 
fiber-optic systems is carried out and the error probability predictions are verified 
experimentally for a four-level PAM system. It is shown that, using multilevel PWM, one 
can achieve very low error rates, provided that a fiber with sufficiently large bandwidth 
is used. Multilevel PAM is more useful where such wideband fibers are not available. 

The mathematical model makes use of the Gram-Charlier series to represent the 
probability density function of the receiver output. Introduced into the Gram-—Charlier 
series model is the effect of intersymbol interference, which is shown to be an 
important factor in limiting the performance of a fiber-optic system. Numerical examples 
are provided to assess the performance of 4—-level PAM and PWM systems. 

The experiment involved the transmission of a four level pseudo-random PAM 
sequence, operating at 4.023 Mb/s, over a 0.5 km. multimode fiber. Error rates of the 
system were measured as a function of the avalanche gain and the optical power incident 
on the photodiode. The experimental values were seen to match closely with the 


theoretical predictions. 
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CHAPTER | 
INTRODUCTION 


TeEIBER = ORMICISYoreEMs 
The potential of optical fiber as a versatile transmission medium for 
communication applications is well recognized [1],[2]. Attractive features of optical fiber 
compared to conventional systems are: 
ie larger bandwidth and lower loss, 
2. reduced size and weight, 
3. higher system channel capacity, 
lower basic material (glass) cost, 
freedom from electromagnetic interference and pickup, 


higher temperature capability and 


eo Ol 


greater security since it is difficuit to tap. 

The attractive features of fiber-optics as a means Of Communication have been 
recognized for a long time. but could only be realized after the development of !ow-loss 
fibers in 1976 {3]. Improvements made in the design of low-loss fibers spurred the 
growth of related components for optical—data links, such as, light sources, modulators, 
multiplexers, couplers, connectors, detectors, demodulators, demultiplexers and 
receivers. A number of field trials, eg, those cited in refrences [4],[5], were conducted 
using fiber optical data links, and these performed sufficiently well that many of these 
links are still in regular service. 

The purpose of a fiber optical communication link is the transmission of 
information from one point to another without introducing unacceptable degradation of 
the signal. This degradation normally is measured in terms of the bit error rate (BER) or 
signal-to-noise ratio (SNR) and signal distortion. Either analog or digital modulation can be 
used in fiber transmission. In a digital communication link, the signal consists of a series 
of discrete pulses which must be distinguished from each other upon reception. In an 
analog communication link, the information is contained in the detailed shape of the 


transmitted waveform. 
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The choice between analog and digital transmission depends on which type will 
provide the required quality at the overall lowest cost. Analog intensity modulation is 
suitable for low bandwidth signals and has the appeal of system simplicity. On the other 
hand, digital modulation is highly immune to noise, provided that the system has a large 
bandwidth. Digital modulation is therefore ideally suited for fiber transmission, a medium 
of large bandwidth. For digital transmission, a SNR of only about 20 db. yields a BER of 
10-°, a figure accepted by industry for high quality systems. Also the frequency division 
multiplexing equipment used for analog transmission is more expensive than the time 
division multiplexing used for digital transmission. The disadvantages of digital 
transmission are that a wide-band transmission medium and complex encoding/decoding 
circuits are required. Fiber-optics, with its inherent wide-band capacity, is ideally suited 
for digital transmission. 

Digital transmission tends to concentrate on the bit-rates used for telephony and 
derivatives of those bit-rates. This has been done to provide standardised 
modulator/demodulator equipment. Thus preferred rates are 1.5, 6.3, 44.7. 274 Mb/sec 
and 1.1 Gb/sec. There is a trend towards using higher modulation rates; however, the 
present characteristics of light emitting diodes (LEDs) limit the achievable bit-rate to 
about 140 Mb/sec for a line length of about 500 m. when using binary modulation. The 
modulation capability using lasers reaches 560 Mb/sec and operation at 1.1 Gb/sec is 
also reported [6]. 

A method used to improve the transmission capacity of the LED or the injection 
laser diode (ILD) employs multilevel signalling. The relation between the information rate 


and the number of levels is given by [7]: 


R = log,(MN/T) (1.1) 


where: 
R = information rate in bauds/sec. 
M = number of levels, assumed to be equally likely 
T = digit interval in seconds. 
Thus, increasing M from 2 to 4 increases the information rate by a factor of two. One 


could also decrease T in order to increase the information rate; this would, however, 
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necessitate a higher bandwidth system. This has been discussed in greater detail in 
Chapter Ill A comparison between the binary and 4-level modulation schemes, with 
respect to the required optical power, for a system using a commercially available LED 
and a PIN-JFET trans—impedance receiver is given in Fig.1.1 [8]. Intersymbol interference, 
i.e. the effect of pulses from adjacent time slots on the pulse under consideration, has 
been neglected. A real-life link would be subject to intersymbol interference and this 
would impose a penalty of a few extra decibels of peak power. The penalty for using a 
larger number of levels is increased intersymbol interference, and thereby greater power 
penalty. Thus, though 8-level modulation would seem to increase the information rate by 
three times, in actual practice the increase is much less due to increased intersymbol 
interference. Another disadvantage of using a larger number of levels is that the 
electronics becomes more complex. Four-level modulation seems to be a good 
compromise between binary and 8-level (or higher) modulation for obtaining increased 
information capability without creating too much intersymbol interference and circuit 
complexity. 

The transmission capacity for an ILD might be doubled, as for the LED, by using 
four-level modulation. This would, however, require a very high or a very well defined 
extinction ratio (on-to-off ratio of light levels). Both these requirements are beyond the 
Capability of currently available ILDs. Also, bit rates of 560 Mb/sec and 1.1 Gb/sec will 
not be readily usable until gigabit logic becomes available. 

Another possible modulation scheme, suitable for multilevel transmission, is the 
Pulse-Width-Modulation (PW/M) scheme. PWM is suitable for high bandwidth systems 
requiring extremely low bit error rates. The SNR of PWM can be improved by decreasing 
the rise and fall times of the optical pulses; i.e, by using a higher bandwidth system. This 
is not possible with PAM systems. The requirement of short rise and fall times makes the 
electronics necessary for PWM schemes rather complex. The finite rise times of 
commercially available LEDs, lasers, fibers, detectors and the associated circuitry limit the 
achievable error probability for a PWM system. Multilevel PWM systems provide 
increased transmission capability compared to binary PWM systems, in the same manner 


as for PAM systems. 
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Figure 1.1 Comparison of binary and four-level modulation schemes at (a)140 Mb/sec 
and (b)70 Mb/sec, using a LED-PIN configuration. The curves show minimum 
peak power required to achieve BER of 10-°, versus fiber length. 

The solid line, ——————, is for a fiber with modal dispersion of 4 ns/km. 


The broken line, - - — -, is for an ideal fiber without modal dispersion. 


an © 70722 


dtenad o3diG 


‘% — — ae fe ad dante @ ra Se 
tual 4 é es 4 ’ at v 


— =< a 7 x —o 
‘ont \ ga & 423 a | 
. * \ | - e~ 
. 
= 74 , - } 
eve ‘_ Sees Lo 
: : =>. ‘ p 2 
i ~ 
anaes awd ae 
7 {ard J SY” mam) Vara — + ee shee, 
r 20 0 


am toa eho eek i argh, Pa Wireapeed ho 

: atl eos eA een weap oe ¢ 
. - ' bd 

mtv pei sligeaee nese a Tie cen in 


: sella etl 
dipemestel hee a aan 


wundine ‘one : 


ales 
- 


_ aan 


1.2 OPTICAL-COMMUNICATIONS LINK 
A generalized diagram of a digital fiber-optic communication link is shown in Fig. 
1.2. This diagram applies both to the binary and multilevel modulation schemes. The 


following sections describe each of the stages shown in Fig. 1.2. 


1.2.1 Signal:Shaper/Encoder 

The electrical signal is first fed into an encoder/signal shaper. In a digital system 
the encoder stage detects the incoming data and regenerates and retimes the symbols 
appropriately for the optical driver. Normally, fiber-optic systems employ intensity 
modulation using various schemes such as; pulse code modulation (PCM), pulse-position 
modulation (PPM) and pulse amplitude modulation (PAM). PPM has the disadvantage of 
requiring a large spectrum bandwidth and a jitter free clock. PAM is attractive for 
band-limited systems, whereas PWM is attractive for systems with wide bandwicth. 

There are also a large number of codes possible, such as: bi-phase (Manchester), 
delay modulation (Miller), binary return to zero (RZ), binary non return to zero (NRZ) and 
2-level alternate mark inversion (AMI) [9]. Choosing a code involves consideration of the 
following: 
ii, Power constraint on ILD: The code should have a small duty cycle so as to avoid 


overdriving the ILD. 


N 


Error Monitoring: The code should have provision to detect and correct errors 

occurring during transmission. 

S Dispersion in the fiber (pulse broadening). The code should be such that pulse 
broadening will not result in too much intersymbol interference. 

4. The circuitry required to generate the code should be easy to design. 

Other important considerations are clock regeneration and baseline wander. 

The signal source normally provides an inherent clock along with the message 
sequence. The clock can be synthesized from the message signal at the receiver end. To 
achieve proper clock regeneration and to facilitate the power output control for laser 
transmitters, it is necessary that the data stream be sufficiently random. Scrambling 
Circuits to achieve randomization of the data scheme and descrambling circuits to 


recover the original data are readily available. 
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Figure 1.2 Block diagram of a typical optical fiber communication system. 
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Baseline wander (i.e., the shift in the dc level of a long pulse sequence), arises in 
digital systems for ac-coupled transmission links using RZ and NRZ codes. When a single 
positive pulse passes through an ac-coupled system the output pulse has a long tail of 
opposite polarity. For a sequence of pulses passing through an ac-coupled system, the 
effect of the tails of individual pulses can accumulate to cause the baseline wander 
effect. Baseline wander can be reduced by using bipolar coding. The price one has to pay 
for this is the transmission of a ternary signal (containing the same information as a 
regular binary signal). A better method to reduce baseline wander is to use a scrambler. 


Fig. 1.3 shows & number of coding schemes commonly used. 


1.2.2 Source Driver and Optical Source 

The encoded signal, is applied to the source driver. The driver modulates the 
current flowing through the optical source to produce the desired optical signal. The 
principal requirements of the source are: 

a. good linearity of the output optical power vs. input current characteristics for 
faithful reproduction of the electrical signal, 

b. narrow spectral bandwidth, 

eC] high optical output at low current density, 

d. small emitting area to couple the radiated power efficiently into the fiber, 

é. high frequency response and 

ue long lifetime even with high current density. 

Along with these requirements is the additional constraint that these sources 
should emit radiation in the spectral region where the fiber attenuation is low - 0.8 to 0.9 
um. and 1.0 to 1.6 um. Diodes that emit stimulated radiation are called laser diodes, and 
those that emit incoherent radiation are denoted as light emitting diodes (LEDs). These 
diodes are made of single-crystal semiconductor materials (GaAlAs and GaAs for the 0.8 
to 0.9 um. speciral region and InGaAs or InGaAsP for the 1.0 to 1.6 um. region) and 
usually employ a heterojunction structure, where the p- and n- type layers of different 
energy bandgaps are combined to produce mode confinement. 

ILDs are characterised by a threshold current, |4, , for stimulated emission as 


shown in Fig. 1.4 [10]. Above |,, the carrier lifetime is very short (thus allowing high 
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Figure 1.3 Various coding schemes employed in fiber-optic communication systems. 
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Figure 1.4 Power emitted from an LED and from an ILD as a function of direct current. 
The knee in the ILD curve is at the threshold current which depends on the 
device and its construction. The nonlinearity in the curves at high currents is 


due to junction heating. 
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modulation rates), the spectral emission narrows down to a few angstrom units of 
wavelength, the beam becomes highly directional and the quantum efficiency is very high. 
The LED, on the other hand, has a broad emission spectrum (typically about 100 A° units), 
a broad-angle beam and a modulation capability limited by the carrier lifetime of 
approximately 1 ns. under optimum conditions. 

Although both CW laser diodes and LEDs can emit several milliwatts of optical 
power, their applications differ. The LED is useful for systems requiring limited 
bandwidth capability and short-to-moderate transmission distances. The CW laser diode 
is useful for high data rate and long-distance applications. Being a device with a 
temperature-—dependent threshold, the laser diode is somewhat more difficult to use than 
the LED. Its reliability and lifetime are also not as high as that of the LED. Keeping all the 
above considerations in mind, one can then select the optical source best suited for a 


particular application. 


1.2.3 Source-Fiber Coupler 

The purpose of the source-fiber coupler is to efficiently launch the optical 
power into the waveguide. Source-to-fiber coupling techniques include direct butt 
coupling, the use of microlenses and of tapered or bulb-ended fiber pigtails. Despite the 
many different coupling techniques available [11], the butt-joint connector is the one 
most commonly used. A butt—joint tolerates lateral and axial displacements in the fiber of 
about 10% of the core diameter and angular misalignments of a few degrees, without a 


drastic increase in coupling loss. 


1.2.4 Optical Cable 


The optical cable transmits the optical signal from the transmitter to the receiver 
either over a single fiber or over a fiber bundle. The principal requirements are low loss 
and low dispersion but other criteria) such as the dimensional characteristics, modal 
interaction, bending radius and economic considerations, may be equally important The 
fiber types that can be employed in wide-band networks are the graded-index (Gl) and 
single-mode fibers. The advantage of the single-mode fiber is its extremely large 


bandwidth of many Ghz km. This is obtained by using a very small core diameter, which 
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imposes extremely tight tolerances on splicing and connector techniques. GI fibers have a 
much larger core diameter than single-mode fibers but still have fairly large bandwidths 
(in excess of 1 Ghz km). GI fibers, because of their ease in splicing and making 
connections, are much more commonly used than are single-mode fibers. 

Depending on the fiber, source and detector characteristics and the total system 


length, it may be necessary to regenerate the optical signal by the use of repeaters. 


1.2.5 Repeater 

The repeater acts as a regenerative system element and is designed to enhance 
the amplitude and the shape of the signal degraded during transmission over the optical 
fiber. It consists of a photodetector, amplifying and reshaping circuitry, optical source 
and coupler. Conceptually, the repeater can be described as a back-to-back 


receiver-transmitter combination. 


1.2.6 Fiber-Detector Coupler 


The purpose of the fiber-—detector coupler is to efficiently couple the optical 
signals radiating from the fiber into the photodetector. It is designed to minimize the 
reflective losses at the fiber-detector interface and to match the respective 


cross-sectional areas. 


1.2.7 Detector 

The detector converts the received optical power into an electrical current. It 
must be able to follow the signal emerging from the fiber both in amplitude and phase. 
These devices should be able to operate in the range of optical frequencies between 0.8 
to 0.9 um and between 1.0 and 1.6 um. In most optical fiber systems the photodetector 
used is a p-i-n or an avalanche photodiode, though in principle a phototransistor or 
photomultiplier could be used. In the 0.8 - O93 um band these photodiodes are made 
from silicon which provides an excellent spectral match to the emission from sources 
made from the AlGaAs system. In the 1.0 - 1.6.um band, detectors using germanium have 
been made but they have higher leakage current and higher excess noise than their silicon 


counterparts. Devices using quaternary compounds such as GaAlSb, GaAlAsSb and 
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InGaAsP are also being developed at these longer wavelengths and they promise 


characteristics similar to those of silicon detectors at shorter wavelengths. 


1.2.8 Amplifier and Signal Shaper-Decoder 

The amplifier enhances the electrical signal generated by the detector and 
increases it to a level at which it can be reshaped for proper further use. The amplifier 
must have low distortion and its bandwidth should be larger than that of the signal. The 
signal shaper and decoder finally converts the electrical signal into a form suitable for 


use. Its design is a function of the intended application. 


Poelhes |S OBJECTIVES 
The major objective of this thesis is to correlate the experimental results 
obtained for a 4-level PAM system with the values theoretically predicted using the 

Gram-Charlier series model. In keeping with this major objective, the following points are 

considered: 

fe To describe and compare various models used for evaluating error probability in 
fiber-optic systems. 

2. To show the suitability of the Gram-Charlier series model for error probability 
analysis of fiber-optic systems. Also to develop this model so that it can be applied 
to multilevel PAM and PWM systems. 

3. To incorporate intersymbol interference into the Gram-Charlier series. 

4. To design and test a 4-level PAM system using an optical fiber as the transmission 
medium. The error probability is evaluated as a function of the average optical 
power reaching the receiver and as a function of the average avalanche gain of the 


APD. 
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1.4 ORGANISATION OF THE THESIS 

The remaining chapters of this thesis have been organised as follows: 

‘Chapter | discusses the different techniques available for estimating the error 
probability for fiber-optic systems. Considered in this chapter are the various noise 
sources existing in a fiber-optic link and their effect on the receiver sensitivity. Special 
emphasis is given to the Gram-Charlier series method because of the computational ease 
with which error probability analysis can be carried out using this method. Incorporated 
into the Gram-Charlier series is the intersymbol interference noise term. 

Chapter Ill describes multilevel PAM systems. The Gram-Charlier series approach 
is used to compute the optimum power levels and thresholds needed for minimum error 
probability. The effect of intersymbol! interference on the error rate of a 4-level PAM 
system is analysed. Graphs of BER vs. the average optical power and BER vs. the mean 
avalanche gain of the detector are also plotted and the results are analysed A 4-level 
PAM system employing 2 fiber-optic link nas been designed, fabricated and tested for a 
transmission rate of 4.0 Mb/sec. A detailed analysis of this experiment is described in 
Chapter Ill. The error probability for different power levels and avalanche gains are 
obtained experimentally and these are compared to the theoretica! values. 

Chapter IV deals with the analysis of multilevel PWM systems. The comparison of 
PWM and PAM systems is also given in this chapter. Multilevel PW/M is shown to be an 
attractive alternative to multilevel PAM, especially for fiber-optic systems. 

Chapter V presents the overall conclusions of this thesis and gives some 


recommendations for further research. 
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CHAPTER II 
ERROR PROBABILITY CALCULATIONS FOR FIBER-OPTIC RECEIVERS 


A key element in any optical fiber communication system is the receiver. The 
basic purpose of the receiver is to convert the modulated light incident upon it into an 
electrical signal, from which the information impressed on the light at the transmitter is 
recovered. A digital receiver essentially consists of a photodetector and preamplifier, an 
amplifier, a filter and a threshold detector. 

The incident light is converted into an electrical current by the photodetector. The 
preamplifier enhances the low-level electrical signal without introducing excessive noise. 
The filter is used to limit the bandwidth (and thereby eliminate the out-of-band noise). 
The threshold detector checks if the output signal is below or above a certain threshold 
and thereby decides whether a ‘ONE’ or a 'ZERO' was transmitted '. Depending upon the 
overall noise of the system, the threshold detector will make a certain number of errors 
in decision. The characterisation of the noise sources and their effect on the bit error 
rate forms the basis of this chapter. 

Several methods exist for the analysis of error probability in optical—fiber 
receivers. The analytical approach of Personick [12][13] provides the most 
comprehensive explanation of the various noise sources involved. A slightly simplified 
approach, suggested by Smith and Garret [14], is described here in detail. Some of the 
other methods, more suitable for computational purposes, are the Chernoff Bound 
method [15], the Characteristic Function method [16],[17], the Monte Carlo simulation 
method [18] and the Gram-Charlier Series method [19]. The relative merits and demerits 
of these methods are discussed here. A detailed analysis of the Gram-—Charlier Series 
method is given because it will be used for the subsequent analysis of multilevel PAM and 
PWM systems. Introduced into the Gram—Charlier Series are the intersymbol interference 


and the thermal noise terms. 


1A binary PCM transmission is assumed here. The two signal levels are denoted by ‘ONE’ 
and 'ZERO’ which correspond to the higher and lower levels, respectively. 
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2.1 NOISE SOURCES IN A FIBER-OPTIC SYSTEM 

Traditional communication systems involve transmission of electrical signals 
through wires and cables or via microwave links. The system in all these cases is subject 
to noise whose mathematical distribution is Gaussian. Such a noise arises from a large 
number of independent small fluctuations in the transmission medium. The Central Limit 
Theorem asserts that a system consisting of a large number of independent noise 
sources has a Gaussian distribution. Because of its simplicity, most of the mathematical 
analysis in traditional communication systems has centered around Gaussian noise. The 
Gaussian model is often used to describe an optical fiber system, but such an approach is 
inaccurate. Lignt detection is based upon photon counting, which is subject to statistical 
fluctuations of a Poisson nature This Poisson noise causes the probability density 
function (PDF) to depart from the symmetry of a Gaussian density function and 
introduces & certain skewness in the PDF. 

The basic photodetector can be modelled as a current source in parallel with a 
Capacitor and a resistor. This is shown in Fig. 2.1, where Cq represents the overall 
Capacitance (including the photodiode, packaging and lead capacitances) and Ry 
represents the biasing resistor connecting the detector to the biasing voltage. Consider a 
light pulse of period T seconds incident on the photodetector. The action of the 
photodetector is to absorb the incident liaht power pit), and to generate hole-electron 
pairs at times {ta} . These hole-electron pairs separate under the influence of the 
internal electric field, thus resulting in a displacement current i. (t). The average value of 


this current is given by: 


<i,(t) > = Rp(t) (2.1) 


where: R = responsitivity of the detector 
= emitted current / incident power . 
The average number of hole-electron pairs, A , that are produced during an 


interval T is given by [ 12}: 


A = (n/hv) fo’ p(t)dt = (n/hv)E ee 


where: 
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ig(t) Cy Rb 


JAYS ene biasing resistor 

Cy —> photodiode capacitance 
+packaging and lead capacitance 
+stray capacitance 


i,(t-—» photodiode current 


Figure 2.1 A simple model of a photodetector along with its associated biasing resistor. 
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n = quantum efficiency of the detector (ie. the fraction of incident photons which 
release hole-electron pairs.) 

E = energy of pulse, pit), over interval T. 
Equations (2.1) and (2.2) provide a simplified expression for the current output and the 
number of hole-electron pairs produced by the photodetector. A complete description 
of these is only possible by considering the noise arising due to the randomness of the 
actual number of hole-electron pairs, N, generated by the detector. In addition, the time 
of occurrence of these events, i } , iS also random. Such a system can be modelled as 
a Poisson process. 

The probability of a photon arriving in a given time interval is proportional to the 
intensity of the incident light. Thus, a strong signal corresponds to a large number of 
similarly sized, but randomly timed events. 

Fig. 2.2 graphically depicts the photon-counting process. The waveform of pit) is 
assumed, for illustrative purposes, to be triangular. Assume that the electron-hole pairs 
are generated at random times ha For a perfect optical detector, the voltage response 
across the load resistor, for a single electron-hole pair, should be an impulse. Fig 2.2(c) 
gives the output of such an optical detector when a single photon has been detected. Fig. 
2.2(d) shows the same output considering the finite response of the detector, hit). The 
overall pulse output, vit), resulting from the responses to all the electron-hole pairs [as 
shown in fig. 2.2(e)] has a smooth mean envelope pulse, which follows the mean arrival 
rate for photons, but has superimposed upon it a large, signal-dependent noise. The 


pulse output, vit), can be described as: 
v(t) = a,h(t-t,) (2.3) 


where a, are the received levels and hit) is the response of the detector. The value of 
a, depends on the number of counts, N, made during a pulse interval. The probability of 


n electron-hole pairs being produced is given by the Poisson distribution: 


P(N=n) = (A"e™) /n! (2.4) 


From Poisson statistics, a mean number, , , of photons is expected and there is a 


random fluctuation of N about that mean from pulse to pulse. This signal dependent 
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(e) The detected optical pulse which is the sum of the responses 
for each electron-hole pair generated. 


Figure 2.2 The Photon-counting process in a photodetector diode. 
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noise, called quantum noise, is shown in Fig. 2.2e. 

In addition to the quantum noise, there is also the noise due to random gain 
mechanisms of the avalanche photodiode. Such a noise is often termed excess noise. 
Despite the excess noise, avalanche multiplication of the electron-hole pairs can increase 
the sensitivity of receivers. Other noise sources existing in fiber optic detectors are the 
dark current and the thermal noise. There are actually two types of dark current noise. 
The first, called simply dark current, consists of hole-electron pairs which are thermally 
liberated in the p-n junction and which are multiplied by the avalanche gain. The other 
dark current noise, referred to as /eakage current, bypasses the drift region and 
experiences no avalanche gain. Fig. 2.3 shows all these noise sources and, in addition, 
other noise sources in a fiber-optic system, such as beat no/se for an LED, device no/se 
of the active element, surface /eakage current and the therma/ noise of the amplifier. 
Beat noise occurs for an incoherent source (such as an LED) due to beats between 
spectral components. Beat noise is typically very small and is normally neglected in error 
probability and sensitivity analysis. The dark current and the surface leakage current can 
be reduced by a proper choice of the detector material and careful design of the 
detector as well as the preamplifier. The quantum noise, the dark-current noise, the 
surface-leakage-current noise and the excess noise are grouped together as shot noise, 
which is characterised by Poisson statistics. The signal-dependent shot noise results in a 


departure of the PDF from symmetry (i.e., causes skewness). 


2.2 THE ANALYTICAL APPROACH TO DETERMINE ERROR PROBABILITY 

The analytical approach of Personick [12],[13], as modified by Smith and Garrret 
[14] is considered here. Fig. 2.4 is an equivalent circuit for the receiver. The receiver 
consists of a photo-detector with biasing circuit, an amplifier, an equaliser and a 
threshold detector. The photo-detector is modelled as a current source, i, (t), in parallel 
with the capacitance, Cy . The biasing resistor is represented by Rp . The amplifier is 
modelled as an ideal high gain, infinite input impedance amplifier with a shunt capacitance 
and Petetance (Ca and Rp ) across the input. The Johnson current noise (thermal noise) 
source due to Ry, is modelled as ib (t). The amplifier current and voltage noise sources 


are i,{t) and e,(t) respectively. All these noise sources are assumed to be white, Gaussian 
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LED or LASER DIODE 


AVALANCHE PHOTODETECTOR 

AMPLIFIER 

Electronic 
Gain 


Photoelectric Avalanche 


7S effect Gain 
Background 
radiation 
noise . 
Beat from (1) Quantum or (i) Excess noise (1) Thermal 
incoherent Poisson fluctuation -due to random -input resistance 
source Roice gain mechanism (2) Device 
(2) Leakage current (2) Dark current -active element 
(3) Surface Leakage 
Sarg nt current 


Dark current noise 


Figure 2.3 Noise mechanisms in a fiber-optic communication system. 
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and uncorrelated. 


Assume that the optical signal, p(t), incident on the detector is a digital pulse 


stream ?, represented by: 


p(t) = 2 b,h,(t-nT) (2.5) 


where: 
T = bit time slot 
h p (t) = pulse shape; 
or = energy in the n th pulse, which can take two values; depending on whether a 
‘ZERO’ or a ‘ONE’ was transmitted. 
The probability of producing exactly N counts during time interval (t,, tp+T) is given 


by (2.4). 
PLN, (to, totT)] = ( aXe“) MN! (2.6) 


where: 


A= fygo™ A(t) dt 
and where X(t) is the average rate of electron-hole pair production. 


NER) = (qyyAswal GO) Se Ne (2.7) 


A, denotes the dark current counts per second. 


The average detector current output, <i . (t)>, is given by: 


<i,(t)> = (nq/(hv)) <g>p(t) + <@>aro (2.8) 


where: 
<g> = average avalanche gain of the APD 
gq = electron charge 
Neglecting the dark current noise term in (2.8), the average voltage at the equaliser 


output is: 


2Miultilevel systems are considered in the next chapter. 
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<Vous(t) > = (na/(hv) ) <g>Ap(t)*hy(t) *heg(t) (2.9) 


where h bt) is the amplifier input circuit current impulse response given by: 
hit) = ZF [(V/R_) + july] (2.10) 


and where: 
heg (t) = impulse response of the equaliser 


The average output voltage can also be described in series form as: 


Vou(t) > 3 i bho (tn!) (27151) 


where hoy (t) represents the output pulse shape. 
Superimposed on this average voltage are the various noise sources described 
earlier. 


Vout) ca <V out (t) > ar n(t) 


= E bahou (t 27) 


(2a) 
where nit) is the noise voltage at the equaliser output. 
The variance of nit), N, is defined as: 
2 
N = root GO is = veo (t)> or eek es (2.13) 
Assuming all the noise sources are mutually independent: 
N = <n®(t)> + <nf(t)> 
+<rf(t) > + <nk(t) > (2.14) 


where 
ng(t) = The output noise voltage due to the shot noise current, (t), 
produced by the detector, 


ie The output Johnson noise voltage due to the resistor Ro ; 
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n itt) = The output noise due to the amplifier input current noise 
source, i ,(t), and 


net) = Output noise voltage due to the amplifier input voltage 


noise source, e, (t). 


2.2.1 Shot Noise 

Personick [12] has evaluated the shot noise term, nt). at the equaliser output as a 
function of time within the bit time slot, T. This results in a comprehensive but rather 
complicated analysis. Smith and Garret [14] make the assumption that the shot noise is 
constant within the bit time slot. The resulting analysis is much simpler and is still quite 


accurate. The shot noise expression is given by [14]: 


<n,*(t)> = 2q<ip>:<g?>ByRFA” (2.15) 
where: 
<i O77 is the mean unity-gain photocurrent over the bit time T. Contributions from 
neighbouring pulses may affect <ig mT if they spread out into the bit time under 
decision. 


By is the noise equivalent bandwidth, defined as: 


2By = (1/R#) [SS |Heg( f) HAC 4) [?df ] 
=(1/R%) |hy/(Ang<g>) |? SZ |Hou(f) /Hp(£ ) [af eight 
where H(f) is the Fourier transform of hit). 

The shot noise expression of (2.15) will now be evaluated for the worst 
intersymbol interference case, i.e., when all neighbouring pulses are ON. The two possible 
cases, i.e. an ON pulse or an OFF pulse at the decision time, will be considered separately. 

When there is an ON pulse at the decision time, the average power, Pavg (t), 


during time slot {-(T/2) 7/2)| due to a pulse centered at t = nT is: 


Payg(t) = (bi/t) S472 hy(t—nT) dt ate 


Therefore the mean unity gain photocurrent, 
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<ip>r,on = 3 (na/hv) (b4/T) f-y'4? h,(t nT) dt 
=(nq/hv) (b,/T)E f-%%h,(t nT) dt 
=(nq/hv)(b,/T) f= h,(t) dt + nq/hv) (b,/T) (2.18) 


When there is an OFF pulse at the decision point: 


<to>n, ort = 2 (na/hv)(b,/T) f= h,(t—nT) dt 
= (nq/hyv) (b,/T) (1-7) (2.19) 


where: 


y = SHR h(t) dt 


Equations (2.18) or (2.19), when substituted into the shot noise expression (2.15), give the 


worst case shot—noise, <n g *(tl>ony OF <n g *t)>orr respectively. 


2.2.2 Thermal Noise 

The thermal noise contribution to the output noise voltage arises from the bias 
resistor, <Vp’>, the amplifier noise current, <V, ?> and the amplifier noise voltage, <\p'>. 
These are expressed as: 


<v2> = (4k0/R,) ByRtA* 
<i> = 25,ByRiA* 


(2.20) 
<V2> = 28,Bi,A" 
where: 
PB es (Heat tmos 
= |(hv) /(Ana<e>) [22 |How (f) /CHp(f ) Fe (f ) ) [Pdf (2.21) 
and By is as defined in (2. 16). 


These results are derived from the well known formula: 
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Sou(f) = Sin(f) |H(f) |? (2.22) 


where: 
Sout (f) = Output power spectral density. 
Sin (f) = Input power spectral density. 
H(f) = filter response. 
The overall thermal noise term is, therefore: 
<V,2> = <V2> + <p> + <Vi> 


= ((2k@/R,) +S; 2ByR3A? + 25,By A? 


(2:23) 
The overall noise term, including the shot noise and thermal noise is: 
<v2> = [((2k@/R,) +S) A? + q<ip>s<g?>A?]- 
2ByR2 + 28,Bi,A? 
= [ (2k0/R,) +Sy+q<ip>rg7"*] - 
A*|(hv) /(Ang<e>)|?S= |Hou(f) /H,(£) |?df 
(2.24) 


+ SpA*|(hv) /(Ana<g>) |?f2 [Hou (£) /CH)(f ) A(t ) ) |?df 
Here, the approximation <g’> = < g>?** has been used [20]. A more accurate expression 
is given later, but to simplify the calculations, it is not used here. 

In order to make the bandwidth integrals of (2.24) independent of bit time, T, and 
a function of only the pulse shape, the following dimensionless time and frequency 
variables are introduced: 
f= t/T and ¢= fT. 


Therefore: 


HL (9) = Hf) and Huw'(¢) = Hu(f)/T (2.25) 


H p | Y ) and Bout (¢ ) depend only on the shape of H p'#) and H out (w ) and not 
upon the time slot width, T. 


Let us define the following integrals: 
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Ie = S2 |How(¢) /Hj(¢) |? de 


I, = 2 ¢* |Hou(¢) /Hp( 9) |? de 
Then: 
JemelHoa ty 7 Hat) idt = Ti, 


and 
—- |How(f) /(H,(f ) H,(f)) |?af 
= J \How(f) /H(f ) |?|RyptjoCy|7df 
= Reis + (2nCy) ShyAt 
Substituting (2.28) and (2.29) into (2.24) gives: 
<V2>= [ (2k /R,) +81q<ip>1<g >?) [hv /( a<g>) |*Tlz 
+S,/hv/(nq<g>) |?] RITIZH 2nCp) *13/T] 
=( hun) *{ <g>*<ip>yTIe/q + Z/<g>*} 
where Z is the thermal noise term given by: 
Z = [St+(2k6/R,) + Sg/Ry) ] Tle/<g>? 
+ (2nC) *Szl3/(T<g>*) 
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(2.26) 


(227) 


(2.28) 


(2.29) 


(2.30) 


(25331) 


From the above equation one can clearly see that increasing R+ and decreasing Cr 


lowers the thermal noise term, Z. This is made use of in the high input impedance 


amplifier approach. Other amplifier designs are discussed in Chapter 3. 


Using the above equations one can proceed to calculate the receiver sensitivity 


and the optimum avalanche gain. 
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2.2.3 Receiver Sensitivity 

Although the shot noise, <n . *(t)>, has a Poisson distribution, the inaccuracy 
resulting from using the Gaussian approximation is not large [12]. The other noise term, 
i.e. the thermal noise, is already Gaussian. Therefore the overall noise term, <v,, (t)>, is 
also Gaussian. The signal amplitude dependence of the shot noise is taken into account by 
the variance of the Gaussian distribution. Pon and Corre are defined as the worst case 
values of <v,y *(tl> for ON and OFF pulses respectively. These are obtained by 


substituting <n 5 Itl>_,) and <n Atl>ore into (2.14), ie. 


Ton = V7 >t] 2 (2) >=205(7) >on 


Corr = VyP>t/*| 7) >=cn,g(T) >oFF (2.32) 


Let m, and m, be the mean signal values for an OFF pulse and an ON pulse respectively, as 


illustrated in Fig. 2.5. Assuming that the threshold, Vi, . is set to give equal area 
probability, Pe , for ON and OFF pulses (ie. the equal area criteria), then: 
[ (m,-Ven) /7ox) = [ (Vino) /Oorr] = Q (2.33) 
where 
Py = (27) °SQ" exp(-x?/2) dx = 0-5 erfc(Q/v) (2.34) 


If m, is zero, then from (2.33) 


[(m,-Vin) /on) = (Ven/orr) = Q 


(2235) 
or m, = Q( Cont+Torr) 
=(Q/<g>) (hv /n) {[ <g>***(n/hv) balptZ] 1 
+ [<g>***(n/hv) deal 1-7) IptZ] 1%} 
(2.36) 


Equation (2.36) is the required sensitivity of the receiver for a given P- (and therefore 
for a given Q value). For a PIN diode detector (with <g> = 1) this equation reduces to the 


simple form: 
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Figure 2.5 Probability density curves for an OFF and an ON pulse assuming Gaussian 
distribution. The mean value of the OFF and ON pulse are m, and m, 
respectively. The corresponding variances of the OFF and ON pulses are Corr 
and ION - The threshold, V4, , is chosen using the equal area criterion (ie. 


the shaded areas A and B are set equal). 
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m, = 2Q(hv/n) vo (2.37) 


2.2.4 Optimum Avalanche Gain 


The optimum avalanche gain is found by differentiating (2.36); ie. dbp, / 3<g> =0 


This results in the following equation from which the optimum gain, g opt: is calculated: 
opt? my = (hv/n) (Z/2I2) ( (2-7) /(1-7) ) K (2.38) 
where: 
k=-1-4{ 14+16( (14x) /x?) ((1-y) /(2-7)2) ] 12 
Corresponding to this the minimum value of m, is: 
{ eex! X 1 
a 1+x CPX er mex (2.39) 
Us amiga Q (hw/n)Z I, L 


where: 


L** = [2( 1-7) /k(2-y)]- 
{[0-5( (2-7) /(1-y) ) k4+1] 1 + [0-5( 2-7) k41] 1/43 2 


2.2.5 Remarks on Analytical Methods 


The analytical approach basically makes use of the theory developed by Personick 
[12]. The probability density function, p(x), can be derived from equations (2.30) and 
(2.31). However, an exact analytical solution to these equations is not possible without 
making a few simplifying assumptions. One of the most commonly used approximations 
is the Gaussian approximation. The major advantage of the Gaussian approximation is that 
one requires only the mean and variance of the output voltage. These are normally easy 
to obtain. The Gaussian approximation, however, does not take into account the extreme 
skewness of the avalanche gain distribution. As a result it tends to overestimate the 


optimal avalanche gain and underestimate the optima! threshold level. An alternative to the 
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Gaussian approximation is the Chernoff bound method [15]. The Chernoff Bound 
replaces the exact solution by a simpler solvable formulation that gives the upper bound 
solution. The Chernoff Bound method is very effective because it gives us the upper 
bound of the error probability. Nevertheless, to estimate how close these bounds are to 


the real solution requires the use of numerical and statistical methods. 


2.3 OTHER APPROACHES TO ERROR PROBABILITY CALCULATIONS 


Consider the fiber-optic communication system depicted in Fig. 2.6. The current 
output of the photodetector is denoted as x(t) and the filter output as y(t). To compute 
the error rate, one needs to know the statistics of the voltage at the decision making 


point, ie. at the output of the filter in Fig. 2.6. The output of the filter is given by: 
y(t) = f2-x(7)h(t—) dr (2.40) 
The error probability can then be computed as: 


Pg = aoPo( y) + aiPily) (2.41) 


where: 
a, and a, are the probabilities of sending ZERO and ONE signals, respectively and 
Y is the threshold. 
“P (y) = P(yIZERO) dy is the cumulative distribution for a ZERO signal 


transmitted. 


= dle P(yt,|ZERO) dy (2.42) 


P (Y) = P(ylONE) dy is the cumulative distribution when 


a ONE signal is transmitted. 


= {2 plyt,|ONE) dy (2.43) 


There are three different approaches to solving (2.40) and (2.41). These are: 
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Figure 2.6 Simplified fiber-optic receiver block diagram. 
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1. Analytical methods 
2. Numerical methods 
3. Statistical methods 
The Analytical approach has already been discussed in great detail. The Numerical 


and Statistical methods are discussed in the subsequent sections. 


2.3.1 Numerical Methods 

Many numerical methods exist for solving (2.40) and (2.41). Normally, x(t) is 
broken up into two componerts, x(t) and Xtp (t), Le. the signal output of the avalanche 
diode and the Gaussian thermal noise generated by the amplifier. If t, is the sampling time, 


chosen so as to maximise yit,), then: 


y(to) = f2 x(t) h(to-7) dr (2.44) 
Since 
x(t) =x,(t) + xa(t) ae 
Therefore 
y(to) a iL X,( T) h(t 9-7) dt 
PS Sel (py igen) ear 
(2.46) 


Vehometeyi hn to) 


The output voltage, yit), is consequently composed of two parts; the signal dependent 


voltage, y. (to) and the Gaussian term, Vin (te 2 The standard deviation of y,, is 
1/2 
On = (4kTBY/R,») (2.47) 


where Bw is the effective noise bandwidth of hit) and Rip is the thermal resistance. 


The effective noise bandwidth is defined as 
B, = (1/|Hol*) So tH(f ) |?df (2.48) 


where H(f) is the Fourier transform of hit) and IH,! is the maximum absolute value of H(f). 


Assuming a rectangular shape for x. (t) and a raised cosine shape for hit), the discrete 


‘The probability distribution of y th’! is Gaussian since x th is a thermally generated 
noise. 
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approximation of the signal y a (t,) is: 


Vara L/h) DL Xqhy-n (2.49) 


where T is the duration of the signal and N is the number of samples used in the 
approximation. The samples xX, are independent and have a distribution dictated by 
equations (2.30) and (2.31) for a time interval 7 = T/N. An approximation for the 
probability density function (PDF), p(x), suitable for computing purposes, was obtained by 
Webb, Mcintyre and Conradi [21]. This approximate PDF for the current at the output of 
the APD is: 
p(x) = {1/[(2m)° of 14 (x-M /or) ]*/7]}- 
exp{ —(x-M) #/[ 20°[ 14{ (x-M) /od) 11} on 
where 
M= (ae - Q/T)G = mean output current 
q = charge of an electron 
T = time interval between pulses 
G= average avalanche gain 
Ne= (nT/hy) Fopt = number of primary electrons 
n = optical conversion efficiency | 
Popt = average optical power in time T 
hv = energy of a photon 
N=V Ne’ Fe (F,- 1) 


o = Gane F. 


e - (q/t)? = variance of the diode current 


F=keff Gt [2 - (1/G)](1 — k ogg ) = excess noise factor 


ef 
K ope effective ionisation ratio 


Since the random variables x, are independent, the density of their sum equals the 


n 


convolution of their respective densities. 
p(ye) = [1/(AT™ hal) 1 {PL x1/( AT: bys) J 
P[ x2/( AT: by-2) J++ °°: + *P[ xy/( AT: ho) ]$ O54) 


The thermal noise is also additive and independent, therefore the overall probability 


density function, ply), of the output voltage, yit), Is: 
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P(y) = plys)*p( yen) (2.52) 


The fastest and simplest way to compute ply) is to transform all partial densities 


using the Fast Fourier Transform (FFT), multiply them and take an inverse FFT: 
p(y) = [1/|(AT) “I(h,) |)- 


Hates | B(x,/A1T* hve) 1} lL pCyenw |D (2.53) 


The Fourier Transform method suffers from numerical errors since it is computed as a 
sum of positive and negative numbers 4 

Direct convo/ution produces less truncation error in ply), because the probability 
densities are positive functions. This method, however, requires an excessive number of 
arithmetic operations. For example, if we have n distributions each having M points, the 


number of operations (multiplication and addition) will be: 


N = n! M’/e2 (2.54) 


For M = 100 and n = 2, N = 10% which is a manageable number. But for M = 100 and 
n= 10, N = 3.5 X 107; which is excessive. Thus the convolution method is of use for only 
limited cases. 

The Gauss Quadrature ru/e method is a numerical technique used to evaluate the 
integrals for the conditional probability densities in (2.41). Most numerical integration 
techniques consist of approximating the integrand by a polynomial! in a region and then 
integrating the polynomial exactly. If the integrand is a complicated function, 
approximating it by a polynomial may be difficult. Often a complicated integrand can be 
factored into a weight function and another function better approximated by a 
polynomial. As an example, if g(t) is a complicated function which has to be integrated for 
t between the limits c and b, it can be factored as wit)f(t) where wit) is the weight 


function and f(t) is a simple function (i.e. one which can easily be approximated by a 


‘Since the regions of interest are the tails of the distribution (error rate of 107%), where 
the various terms consist of small positive and negative numbers, the numerical errors as 
a result of addition will be large. 
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polynomial). Thus, 
fe g(t)dt=foe w(t)f(t)dt=B, wi (ty) (2.55) 


The quadrature rule {Wj t,}, corresponding to the weight function wit), is available 
sometimes in tabulated form *. Often, however, these rules are not available in tabulated 
form. Golub and Welsch [23] have presented algorithms for generating the Gaussian 
Quadrature Rule, provided the moments of the weight function are known or can be 
calculated. This is precisely the situation existing for the fiber-optic system. The Gauss 


Quadrature rules can be used to evaluate the conditional probabilities in equation (2.4 1). 


Poy) = Sy” PC yo[ZERO) dy = SZ f(y) p(y1,|ZERO) dy 


=z, £0 () wi 


(2.56) 
Similarly, 
Py) =f Cw Boy 
where 
1 for y2y 
bal Cy) (2.58) 
0 otherwise 
and 


KOSI S HOGS (2.59) 


The abscissas y; and the weights w, may be calculated by methods of functional 
analysis from the conditional moments Ef(y(t,))Ibo]. Hauk et. a/. [24] have evaluated these 
integrals and thereby derived the error probability. 

Other numerical techniques include the Gram-Charlier Series method [19] and the 
Characteristic Function method [16],[17]. The Gram-Char/ier Series method involves the 
calculation of the cumulants of individual noise sources which can then be added up to 


SThe book "Handbook of Mathematical Functions” by Abramowitz and Stegun provides an 
excellent tabulation of some useful weight functions [22]. 
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compute the overall cumulant, and thereby the probability density function. The major 
advantage of the Gram-Charlier Series method is that it requires very little computing 
time as compared to other methods. A detailed analysis of the Gram-Charlier Series is 
provided in the next section. 

The Characteristic Function method is essentially similar to the Gram-—Charlier 
Series method: instead of adding the cumulants of the individual noise sources, the 
characteristic functions of the individual noise sources are multiplied to obtain the overall 
PDF. The PDF of all the noise sources is described by a trigonometric series. The total 
characteristic function is related to the Fourier coefficients of the periodic probability 
density function. This results in an expression which can then be used to determine the 


error probability for a given bit rate as a function of the transmitted power. 


2.3.2 Statistical Methods 

The only statistical method that will be considered here is the Monte Car/o 
simulation method. The Monte Carlo simulation method is normally employed for 
systems whose parameters have known statistical distributions. The Monte-Carlo method 
starts with the use of the statistical distribution of the system together with random 
number generators to produce a large number of sample systems. These systems are 
then analysed on the computer and empirical results about system performance are 
obtained. There are, however, limitations on the straight-forward application of the 
Monte-Carlo method for fiber-optic systems. Since error rates are typically of the 
order 10-§ to 10°%, at least 10 to 100 samples are needed in this region of low 
probability. This entails the need for as many as 10% samples, which is impracticable. A 
modification known as the importance sampling, solves this difficulty by modifying the 
distribution so that more samples can be taken from the important region. This reduces 
the number of simulation runs by 5 or 6 orders of magnitude. 

A detailed analysis of Monte Carlo simulation using Importance Sampling is given 


in [18]. 
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2.4 DETAILED ANALYSIS OF THE GRAM-CHARLIER SERIES METHOD 


A Gram-Charlier Series expansion for the probability distribution of noise in a 
fiber-optic receiver is considered here. This representation is more general than the 
previously published work of Rice [25] and Personick [ 12]. 

1. It takes into account the randomness of the avalanche gain as well as the amplifier's 
thermal noise and permits the consideration of any generalised processing filters. 
Previously reported calculation methods [15],116],[18] permit thermal noise to be 
included only at the expense of additional approximations and require complicated 
numerical techniques to solve for error probabilities. 

2. Each of the noise sources can be considered separately and, provided they are 
independent, the cumulants of the noise sources can be added to obtain the overall 
cumulant. Thus the effect of intersymbol interference, dark current noise and modal 
noise can be easily incorporated. 

3. The probability distribution of physical noise sources can often be obtained by 
finding the characteristic function and then taking the Fourier transform. There are, 
however, situations where the characteristic function cannot be obtained in closed 
form. In such cases, it is often easier to describe the PDF as a Gram—Charlier series 
expansion. 

4. This method is suitable for application to multilevel PAM & PWM See 

The Gram-Charlier Series involves the expansion of an ‘unknown’ probability 
function in terms of the derivatives of a ‘known’ density function. For a system with a 
large number of independent noise sources, the output is nearly Gaussian by virtue of the 
Central Limit Theorem. The Gram-Charlier Series derives an expression for the PDF 
which exhibits the limiting Gaussian characteristics. For the Gaussian function, only the 
first few terms would exist. As the function becomes more non-Gaussian in shape, the 
number of terms increases. 

If we expand the logarithm of the characteristic function W(s), in a Taylor series, 
the coefficients of this series are, by definition, the cumulants of the stochastic 


process [26]. 


= i n 
In ¥(s) = 2, ArS 7a! (2.60) 
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Central to the Gram-—Charlier Series approach is the use of the cumulants, A, - to 
calculate the PDF due to a combination of noise sources. Thermal noise processes, which 
have a Gaussian distribution, have cumulants only up to an order of two. There are, 
however, processes such as shot noise and impulse noise which have cumulants of 
higher order. These higher order cumulants can be used to describe the non-Gaussian 
character of the process. 

The description of probability density by means of cumulants arises in a natural 
manner if we approximate the functions by means of a series of Hermite Polynomials. 
The Hermite polynomial is generated as the Gaussian function is successively 
differentiated. These functions are orthogonal over the extended real line with respect to 
the Gaussian weighting function. Such a series, known as Gram-Charlier Series, is given 


by {1S} 


p(x) ==) ang, (x) (2.61) 
where: 
¢n(x) = d®/dx™ [(1/(2m)°°) exp( —x® 2) ] (2.62) 
and 
iw rb Neate erat (2.63) 


where U A(x) is the n'th central moment of x. 

Another common practice, convenient for computational purposes, is the 
normalisation of the random variable such that it has a zero mean and unit variance. If m is 
the mean of the random variable x and e& is its variance, and if we denote x as the 


normalised variable, the relation between x and X is given by, 


x = (x-m)/o (2.64) 


The corresponding normalised cumulants are denoted by 


~ 


Onn ne Ol Oke aay (2.65) 


The normalised characteristic function, Wis), is similar to (2.60) with h, replaced by 


geviertan: a Sale? rela: ont deb ataseriag, Tani 
os taeeens? iAey poll by Sakiefot od att ce nae ra aie 
Avri@ cl sowed sshgnt- mee) ae Awol eee 62k rare t : 


res (4) fe a 


‘4 


[Bi fing ils 1 ee ~ ha! 


in a) = os 


(fea) | 
“ igvtarnatn mores tre ot ating 
wit, 2) edeq ing, lanopetugmes, ” tof egapeasiaey “german rome. ‘eefporay ass 
1) cr 4) gonGi-gy Wie he Cbsty oad Seer “ “spminte eb St 1 
arttay 7 Send aw Apert casnelay gn of > rv ap =O) we “ores 
vd navn i pdt: . 
245 
vq bovarah 


(a wegaletah yh dy 08 Shak owirwe bait 


- 
_ 


= & See ae 


40 


2] 


¥(s) = exp{Z A,S"/n! } (2.66) 
The normalised characteristic function, Fs) is also directly given by 
¥(S) = f_2 exp(Sx) p(x) dx (2.67) 


Subsituting the expression for p(x) from (2.61) into (2.67)gives: 
¥(S)=E agf= exp(Sx) d"/dx*[ (1/(2n)°*) exp(-x?/2) ] dx 
oe!) cA oeeXD( Dee) (2.68) 


Substituting relation (2.63) into (2.68): 
Bi re “4 és 
¥Cs)) = [2 {,5"/n! ] exp(S*/2) (2.69) 


Using (2.66) and (2.68) and noting that X, = 0 and = 1 it follows that: 


[EZ maS*/n! Jexp(S*/2) = exp{E, X,S°/n! } 
= exp{i, A,S"/n! } exp(S*/2) es) 
Equation (2.70) can be simplified by dividing both sides by exp(S?/2) and expanding the 
right hand side in a Taylor series. This results in: 
oe jad = exp{E A,S7/n! } 
=1+(£, A,S*/nt) +-0-5] B ApS*/nt ]? +-.- (2.71) 


ow 


From (2.71) the recursive relation between the central moments and cumulants , A, can 


be derived and is seen to be: 


Hy = 1 
Hg = As 


Mo+i => Att + ae BGs AmitMn—m n-3 (2.72) 


Using equation (2.72) and (2.63), one can calculate the Gram-Charlier coefficients ap. 


There still, however, remains the task of calculating the various cumulants. 
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2.4.1 Derivation of Cumulants 
The method used in the derivation of the cumulants is based on the work of Rice 

[26]. Rice obtained the characteristic function of the noise generated in an electronic 

vacuum diode. In order to obtain the probability distribution he required the inverse 

Fourier transform, which was impossible to obtain in closed form. He effected a power 

series expansion in which he could obtain the inverse Fourier transform of the individual 

terms. The resultant series is recognisable as a Gram-Charlier series. The Poisson noise 
or the shot noise is well-suited for a cumulant description of the sort used by Rice. 
The distinguishing characteristics of this type of noise are: 

i It is the result of the superposition of events occurring at random times. For 
example, let us consider a vacuum tube. Assume that the arrival of an electron at 
time t = 0, produces an effect Fit) at some point in the output circuit. If the output 
circuit is such that the effects of the various electrons add linearly, the total effect 


at time t, due to all the electrons, is 
(t) =< aF(t-t,) (2.73) 


where the k'th electron arrives at the time t,& causes effect F(t-t , ) at time t 
(Refer to Fig. 2.7} 
Consider the case of electron multiplication. Instead of only one electron 

arriving att, a k electrons arrive at each instant. The same idea could be used in a 
fiber-optic receiver model. For a PIN diode a, = 1 and for an APD, a re She 
where om is the number of electron-hole pairs produced by a photon arriving at 
time t = tk: 

2 The occurrence times t , are assumed to form a Poisson sequence with an average 
rate \ per second. This condition holds true for the fiber-optic receiver. 

Rice has shown that, for a system having a noise process which satisfies the above 


condition, the n'th cumulant can be given by: 
A, = varf=[F(t)]? dt (2.74) 


A formula similar to (2.74) will be derived for the fiber-optic receiver shown in Fig. 2.6. 
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Current 


F(t) F(t-t,) F(t-t,)  Flt-t,) 


ne ane eee SE 
Olt tenis 


(a) Response to an electron arriving at times t,, t,, t, and t,, 
each considered separately. The arrival of an electron at time t = t, 


produces effect F(t-ty. 


Gain multiplication 


a 


(b) Random gain multiplications. 


The total effect, at time t, due to all the electrons is: 


Do a, F(t ty) 
k=6 


Figure 2.7 Graphical description of the shot noise analysed by Rice. 
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The receiver in Fig. 2.6 consists of an APD, a voltage amplifier and a linear, 
time-invariant filter. The APD and the amplifier are assumed to have an infinite bandwidth 
and hit) represents the total impulse response. Each electron generated optically at the 
input of the APD gives an impulsive current of height (qG) at its output,where q is the 
electron charge and G is the instantaneous random gain of the diode. After amplification, 
an impulsive voltage of height (qGRA) is produced at the input of the signal processor 
which in turn gives an output of [qGRAhit- 7 )] at time t, where 7 is the instant of the 


release of the initial electron. Assuming both P(t) and hit) are zero for t < 0, we can write 


V(T) = $ (qRAG) h( Tr) (2.75) 


where k is the random number of electrons released in the interval [O, T] and 7, Is the 
time instant at which the i'th electron was released. Since the electrons are released 
independently, the characteristic function of each of these electrons can be multiplied to 
obtain the overall characteristic function of Vit) [26]. The characteristic function of the 
ith term of Vit) is 

¥(s) = Efe (sae) H(t 8} 


ay es e (GRAS) R(T) 8p G)- P( 74) dry (2.76) 


where Po () and P| () are the probability distributions of G and T respectively * Also the 


probability distribution of P_ ( t, } is given by 
P(r) = P(1)/fo? P(t) dt Me 
where P(.) is the input power. 
¥( s) = ie a el aRAb(T—7) SG], 
Po(G)[P( 7) /fo P(t) dt] dry 
A So" P(t) dt] fo" M[ qRAh(T-7) S]- aoe 
Po( 7) dt; 


where MI] is the moment generating function of G as defined in [27]. 


6\Note that G is a discrete random variable. 
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The overall characteristic function of k electrons is the characteristic function of 


each electron multiplied k times, ie. 
¥,(S) = (4(8))* (2.79) 
The probability of obtaining k electrons in a time period, T, is 
Ss e"Kreyk (2.80) 


where k is the mean number of electrons transmitted in the time interval [0, T] 


k = (n/hy) fo" P(r) at (2.81) 


and where h is the quantum efficiency, h is Planck’s constant and v is the frequency of 
light. The characteristic function of the output voltage V(T) is therefore given by 
oe tr ' k 
Yr (S) = £, [e#(k)*/e: ](#(8)) 
iva it ¥/k1 
2 es (k ¥(S) ) */k! 


— a @xh5) 
ee (2.82) 
= @(n/bv) fo" P(r) [4 gRAb( Tr) 8} -1] dr 
Since: 
M(S) = %_, M‘")(0)S*/n! 
=1+£  <G>s"/n! ees 
Therefore: 
¥yr(S) = exp{(n/hv)E <G™>fo™ P(r): 
(2.84) 
[ qRAh( T-7) ]7(S"/n! ) dr} 
or 
loge¥yir(S) = » { (7/hv) ( qRA) ®<G®>- 
(2.85) 


So" P( 7) b"( T-7) dr}S"/n! 


The cumulants of the random variable V(T) are: 


7 i 7 a v as ' - , ‘ 
_ 1 a 7 ha 7 
es). : 4 ia) = (ie | _ 
. ed 
7 air peragsIneni shore Vee Ta nt pagan ; 
| | gy 7 - 
“na ) ape 


= 7) a 
(T ONlavisi era att ny basiwmaney BNO voeS eo amen. aay ara at ie 


(fa) ee (OE Th Matt) = 4 


So yoreu eat a av big iraenaod 6 orG aL enol the opeteup at hte 
a eett 21 (TUV Speliny Tei @8 a tar orn) Tcl eeper BAT INH 


ya regs" é f 
ae 2 ‘ D —— a, 
ZW f tay" fry* | * awe 
_ 
he a . | 
P - a 
‘We, + = i 
($$) 5 
eh [t= 18 fen F ead wh Stes 7 aivuw 
- ; 
fac Pak. w= (7 
| * ay 
a U ie) OF - : 
in "E> = a 
. 


i ws o> Sion ieee - (2)ew? 7 
rb tahee | (2 el ; 


(peas) 


mE ay) (oat ee) 
— r te feb | + Kt ¥ 
aes 


45 


An = (n/hv) (GRA) *<G*>- 
for P( 1) bh T-7) dr (2.86) 


Knowing the cumulants, it is a straight forward task to calculate the error probability 
using equations (2.6 1), (2.63), (2.72), and (2.86). Though thermal noise, dark noise, modal 
noise and intersymbol interference are not directly considered, they are taken into 
account by the term Pit) in (2.86). P(t) can be broken up into signal power, dark current 


power, intersymbol interference and modal noise power, i.e. 


P(t) = X(t) + X,(t) + X(t) + X(t) (2.87) 


where 

X q(t) = dark current 

X g (t) = shot noise term = dygit) 

X (t) = Lc eae: interference term 

= 9D, Chaleatine seas 

Xy (t) = aoe noise term 

g_ (t) = signal waveshape at the filter output 
Zauetione (2.74) and (2.75) can be put together to get the overall cumulant 

An = (n/hv) (qRA)*<G*>- 
fot (XatXstXi+X,) h2( T-7) dt 


(2.88) 
= Ana + Ans + Am + Aon 
where 
Note dark noise cumulant 
n e 
Ana = (7/hv) ( qRA) "<G"> oer 
fo Xgh"( 1-7) dt 
Ans= shot noise cumulant 
= »<G*>: 
Ans = (n/hv) ( qRA) sae 


fo’ Xsh"( T-7) dt 


AM intersymbol interference cumulant 
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An -= (/hv) (qRA) *<G*>- 


ip X,h"( T-7) dr (2.9 1) 


AnN = Modal noise cumulant 

Aan = (7/bv) (GRA) *<GP>- 
lee Xyh"( T—1) dr (2.92) 
The overall cumulant due to a number of independent noise sources is the sum of 
individual Cumulants. Correspondingly, the characteristic functions of individual noise 
sources can be multiplied to obtain the overall characteristic function [26]. If we denote 


Vy to be the overall characteristic function due to the shot noise, dark current, 


intersymbol interference, and modal noise terms, then: 
¥x(S). = Yxa(S) Yas S) Fu S) Yor 8) (2.93) 


We make use of this in incorporating the thermal noise and intersymbol interference 
terms into the cumulants. 
Thermal Noise 

Let W4,,(s) be the characteristic function for the thermal noise term. Since 


thermal noise is a Gaussian process, one can expand Vth (s) as 
Fil S) = exp( Tino /e) (2.94) 


Where oon is the variance of the thermal noise. Since the thermal noise is 
independent of the other noise processes, its characteristic function can be multiplied by 


Y x (s) to obtain the complete characteristic function Y (s). 


¥xin(S). = ¥(S) %1(6) 
= exp[ SL A,S*/n! ] exp(o.§5"/2) 
= exp[A,S-+ AzgS2/2) HAsS*/3l ) + + + + 04857/2) J 
= exp[A\S-+( Agta?) S*7/2+(As8°/3! ) + -- | (2.95) 


X-th 


The thermal noise term is seen to affect only the second cumulant. Thus, by 
changing A, to A,+ one , the thermal noise effect can be included, without needing to 


alter the error probability calculations. 
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ntersymbol Interfer 


The cumulant for the intersymbol interference can be derived from the 
characteristic function of the output voltage. If we denote x(t) as the current output of 
the diode and y(t) the voltage at the decision point (refer to Fig. 2.5) then: 


x(t) = Xe (t) + x,, (t) 


th 


Xs(t) ==, exd(t+y) (2.96) 


=signal term which is a doubly stochastic Poisson process. 
here 5 is the time of release of the k'th electron and N is the random number of 
electrons emitted during the time slot. 


The output of the filter, y(t), has a signal component 
Vet) meee ssc) a(t 7) dz: (2.97) 
This has a characteristic function given by 


@y.(jv) = exp} f= €(7)[@z(jvh(t -7) ) -1] dr} (2.98) 


where 


ECT) = [ag(t) /dio] tua(t) HE p4(t) /dio] 


(2299) 


and where x, (t) is the number of electrons at time t and 6 (iv) is the characteristic 


function of the random avalanche gain mechanism. 


@.(jv) = (1-jv<G>)* forse<G> el 0 (2.100) 


The proofs of equations (2.98) and (2.100) are given in appendix 2. One can directly 


derive the cumulants from the characteristic function by noting that: 


Ap = j? d?/dv™[log.( @y,( iV) ) iva 
= j? d®/dv"[ f= €(7) { @[ jvh(t -7) -1] $ dr (2.101) 


Using eas. (2.100) and (2.101) one can readily derive the cumulant to be: 
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An = n! <G>"f_t é(7)h™t-1r) dr (2.102) 


This result appears to be quite similar to the cumulants derived using the Gram—Charlier 
Series in (2.86). If one takes note of the fact that E(t) = P(t)in/hv) and <G™ is closely 
approximated by ni<G> for G > 10; then the two results become identical. 


The intersymbol interference term is just a part of this term. 


NepaeGe jee) hot 7) d7, (2.103) 


where 
ae) = (oa 14( T) /dio] (2.104) 


The equations given above may not only be used to evaluate the intersymbol 
interference separately, bur also to show the similarity between the Characteristic 


Function approach and the cumulant approach. 


2.4.2 Comments on the Gram-Charlier Series 

For computational efficiency it is necessary that the series in (2.60) be truncated 
to a few terms. The number of terms in the series approximation necessary for a 
particular desired accuracy depends on X, which is defined as the number of standard 
deviations from the average, ie. X = I(V-V)/ oO; |. Error rate computations require 
calculations at large values of X (ie. at the tail of the PDF curves). To maintain accuracy at 
large values of X. the truncated series must contain more terms. As an example, the 
truncated series i, ayy’ (V) has less than 1% error up to X = 8, for an average 
laser power of 10 uW. and assuming an integrate-and-dump filter which passes a 
pulse of width 10°' sec. 

The convergence of the Gram-Charlier Series is not immediately obvious but has 
been proved [27] The convergence of the series is not monotonic, but is in groups of 


three terms. Furthermore, the terms a,y'”)(V) do not decrease in importance until 
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nN > No, where n, is an integer which depends on several factors 7. The addition of thermal 
noise causes the standard deviation to change from 9s to (o3+0;2) 1/2 (where a, is 
the shot noise portion and Gi is the thermal noise portion). As a result, the value of X 
changes from (V-V)/ag to (V-V) /(o8+0,2)' thus the n'th term in the 
Gram-Charlier Series decreases by a factor of 1/(1+0,2/08) 1/2) due to the thermal 
noise. This means that the contribution of the higher-order terms (of the Gram-—Charlier 


Series) decreases as the thermal noise increases. 


2.4.3 Error Probability Calculations 

To evaluate the error probability for a fiber optic system it is necessary to 
compute the conditional error probabilities (ie. receiving a zero when a one was 
transmitted and vice versa). Fig. 2.8 shows the probability density function for the 
normalised variable X as well as the non-normalised variable X'. These two variables are 


related by: 


X = (X4m) /o (2.105) 


where m is the mean and o- is the variance. The conditional! probabilities are given by: 
P[E|O] = Suz Pol=') dx’ =o, f Po(x) dx 
i, (n) 
=i Lay [Zz Bn? (x) ] dx 
= 0- §[ 1-erf ( (Ven) /\200) ] 


| (2.106) 
ae Bn Pas ( Vin—mMo) /Ne0) 
Similarly 
P[E|1] = 0 8[ 14erf((Vin-m,) /¥204) ] sone 
+X ag Gn-1( (Veni) /~e0}) 

n=3 i 

where 
2 
Erf (Vin) = (1/( 27) 1/2) fovth e* dx (2.108) 


These two areas must be equated to obtain the near minimum error probability and then 
solved for Vi, . Appendix 3 gives a program which calculates the error probabilities 


and the near-optimum thresholds. 


"For a fiber-optic system, n, depends upon the optical power input, the bit-rate, the 
mean avalanche gain and the thermal noise. 
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PDF curves 


PDF curves 


Vtha 0 Vth,o 


Figure 2.8 Probability density curves for binary transmission (a) using the non-normalised 


variable X’ and (b) using the normalised variable X. 


CHAPTER III 
MULTILEVEL PULSE AMPLITUDE MODULATION 
The main advantage of using multilevel signals instead of binary signals is that, for 
a transmission channel of a given bandwidth, more information can be transmitted. The 


relation between the information rate and the number of signal levels is given by [28]: 


fh = Te log, (M) (3.1) 


where: 
f ph = equivalent bit rate in Mbits/sec. 
M = Number of levels (assumed to be equally likely) 
f s = modulation rate of the fiber in Mbaud. 

Equation (3.1) can be interpreted in two ways: 

I, If f, is kept constant, fs decreases as M increases. This means that a multilevel 
system does not require as high a bandwidth system as is required by a binary 
system., 

2. If fg is kept constant, then f, increases by a factor log,(M). Thus one can increase 
the information rate of the system without increasing the frequency response of 
the system. 

At high bit rates, pulse broadening, which arises due to the finite bandwidth of the 
system, results in intersymbol interference. One can conceivably reduce the intersymbol 
interference by the use of equalisers at the receiver [29] There exists, however, a 
certain bound on the transmission rate which cannot be exceeded despite the use of 
equalisers [30]. To exceed such bounds on the information rate, one is forced to use 
multilevel signals. 

Equation (3.1) seems to suggest that one can increase the information rate to any 
desired value by increasing M. One must realise, however, that (3.1) does not take into 
account the effect of intersymbol interference. Thus, increasing M from 2 to 4 increases 
the information rate by a factor of less than two, for a given error rate, due to the 
effect of intersymbol interference. As the number of levels increases beyond four, 
intersymbol interference becomes even more pronounced. Another important factor is 


that the circuit complexity increases as the number of levels increases. Keeping these 
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considerations in mind, 4—level PAM is a reasonable compromise between increased 
information rate and increased intersymbol interference and circuit complexity. 

In this chapter, the Gram Charlier series will be used to approximate the 
probability density function of the output voltage and to calculate the optimal pulse 
amplitude levels and thresholds for a 4-level PAM system. The analysis used can be 
extended to any general multilevel scheme. Also presented in this chapter are the results 
obtained from an experimental 4-level PAM system that was designed to operate at 
se 8 Mb/sec, using 600m. of graded index fiber, an LED as the light source and an APD 
as the light detector. The experimental results are compared to the values obtained by 
theoretical calculations. 

The 4-level PAM scheme is shown schematically in Fig. 3.1. The message source 
generates a sequence of Statistically independent random variables a; with values a,, a,, 


J 
a; and a,. 


x(t) =£_ a, 6(t-iT) (3.2) 


where T is the time slot available for each pulse. 
The encoder converts the source signal into a sequence of pulses suitable for 


transmission: 
x(t) =£_ a, r(t-iT) (3.3) 


where rit) is an ideal rectangular pulse. 

The encoded pulses, x, (t), suitably amplified, drive the optical source (LED or 
ILD). The optical pulses generated by the source will depend on the bias current, the 
modulating signal current and on the source response. The output of the optical source is 


denoted as: 
x(t) =~. a; ri(t-) (3.4) 


where r,(t) represents the pulse shape at the output of the optical source. 
The optical signal, while passing through the fiber channel, suffers attenuation and 


distortion. The optical signal incident on the photodetector is assumed to be a 4-level 


sequence, given by: 
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Figure 3.1 Block diagram of a multilevel PAM fiber optic communication system. 
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p(t) =~ Pay a(t) (35) 


where 
p(t) is the received optical power 
g(t) is the normalised pulse shape and 
Pmij is the maximum amplitude of the i’th pulse which can take on 4 values; P ,, P 3, 
Fr wand P,. 
The receiver section consists of an APD, an amplifier, a filter and a decision circuit as 


shown in Fig. 3.1. The average output voltage at the filter output, at time T, is [19]: 

Volt) = (n/hv) -(qRA)-<G>-fo" P(r) h(t-7) dr (3.6) 
The n'th cumulant, A , of the output voltage, is given by (2.86): 

An = (7/bv) -(qRA)®- <G®>-fo" P(r) h™(t—7) dr (3.7) 


Equations (3.6) and (3.7) give a complete description of the statistics of the output 


voltage. The output voltage, can then be represented by the Gram-—Charlier series as 


p(x) = Se An?n( x) (3.8) 


where a,, and ¢, are as defined in (2.62) and (2.63). Fig. 3.2 shows the probability density 
function of the output voltage for a 4-level PAM system. Here m,, m,, m,; and m, are the 
mean values of the PDFs, corresponding to the four levels transmitted. X,", X," and X,’ 
are the thresholds chosen to distinguish between the various levels. The areas A,, A,, As, 


A,, A; and A, are the tails of these curves. The error probability is therefore given by: 


P[e] = p(1) {pl2l1]} + p(2) fplile] + plalels 
+ p(3) {p[2|3] + p[4/3]} + p(4)ipl3i4]} (3.9) 


where: 
p(i) = probability of the i'th level being transmitted; i=1,2,3,4 
plilj] = conditional probability that the received signal is interpreted as the 


i'th level, when actually the jth level was sent i,j [1,2,3,4] 
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m,, M., m; and m, are the mean Signal levels. 
X,", KX," and X," are the decision thresholds. 


Figure 3.2 The probability density function of the output voltage for a 4—level PAM 


system. 
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Equation (3.9) can be simplified by assuming that all the signals are equally probable; ie., 
p(1) = p(2) = p(3) = pl(4) = p = 1/4. The equal area criterion is also used (i.e. A, = A, = A; = 
A, = As = A, = A) since it simplifies the calculations and results in error probablities close 


to the minimum [31]. With these assumptions,the error probability reduces to: 

Ple] = p(Azg) + p(ArtAy) + p(AstAg) + p(As) = (3/2)A (3.10) 
Area A, is given by the integral, 

Noyce pS P,(x”) dx" (3.11) 


The other areas can be represented by similar equations. Equation (3.11) can be 


converted to the normalised case (i€. zero mean and unit variance) by using the normalised 


variable 
x = (x"-m,)/o, for P,(x) 
x= (x"-—m,) /oc, for. Po{x) (3.12) 


Using (3.8) and (3.12), equation (3.11) reduces to the following form: 
A, = 0-5[ 1-£rf((x{-m)) /Ne0;) I ser D3 BnYn-1( (xf —m,) /o;) Cake) 


The other areas are derived in a similar manner. 

The derivation of error probability for a 4-level system is a direct extension of 
the derivation used for a 2-level system. Using (3.9) each of the areas can be equated to 
two-thirds of the overall error probability. The optimum signal levels and thresholds 
required to achieve a particular error probability can then be determined by solving (3.13) 
iteratively. However, one can simplify calculations by assuming the PDF's to have a 
Gaussian distribution. For such a case, the Gram-Charlier coefficients, a, , for n > 2 are 


zero. The resulting equations are given below: 


Erf((x{-m,) /éo,) = 1 —- (4/3)P[e] (3.14) 


Erf ((mp-x/) /Méag) = 1 — (4/3)Ple] (3.15) 
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Erf ((x—m,) /eo2) 


tara KEZS)E eh (3.16) 


Erf ( (mg 5-x) /Veos3) 


1 — (4/3)Ple] 


(3.17) 
Erf ((x§-ms) /Veos) = 1 — (4/8)P[e] (3.18) 
Erf((m,-xg) /véo,) = 1 — (4/8)P[e] 

(3.19) 


The o 's in equations (3.14) to (3.19) represent the variance of the Gaussian 


PDF's. From Chapter 2, this is seen to be: 


a, = bh + Ag (3.20) 


where 
o,,, = thermal noise voltage of the system. 


Ai = second cumulant for the ith signal level, having units of volt?. 


3.0.1 Intersymbol Interference 
The intersymbol interference can be incorporated into the Gram-—Charlier 


coefficient for the 4-level case, in the same manner as it was for the 2-level case in the 
previous chapter. To do this, one needs to change the expression for the power input, 
pit), in (3.6) and (3.7), to include the intersymbol! interference terms. Making the 


assumption that only the pulses in the adjacent time slots create intersymbol interference. 


pit) is seen to be: 


‘The intersymbol interference due to pulses in time slots not immediately adjacent are 
very small and can generally be neglected. 
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p(t) = P.g(t-T) + Pog(t) + Pye(t+T) (3.21) 


where P,git) is the signal term and P_,g(t-T) and P.,g(t+T) are the intersymbol interference 
terms from the adjacent time periods. The function, g(t), represents the pulse shape. P_, 
represents the maximum power level of the pulse just before the time period in 
consideration and P,, represents the maximum power level of the pulse just after the 
time period in consideration. P, represents the maximum power of the pulse during the 
time period under consideration ie. (-T/2,T/2). It is assumed here that the maximum 
power for each pulse occurs at the same position in the time slot. 

Expression (3.21) can be substituted in (3.6) and (3.7) to get the overall noise, 
including the intersymbol interference term. There are actually 64 (43) possible 
combinations of levels for intersymbol interference for the 4-level case. One normally 
considers the worst case intersymbol interference ie. lowest signal level flanked on both 
sides by the highest signal level and vice versa. The four power levels will be denoted by 
Pairs. and Pe 

Thus, when the lowest signal level occurs between the two highest levels, the 


expression for the n'th cumulant, Aj, is: 


a a lenVaho.) CGR ANe maaan 


cg ARQ Grea EEA I ECAP AUS BAKE ares, 


Assuming that g(t) is a Gaussian pulse of unit height and standard deviation, o , equation 


(3.22) can be rewritien as: 
An = (n/hv) (qRA)® <G?> 
SEB [PY ett 720" 4.4747) “/eo® 


Sh (323) 
+P,e7 /2¢ ]- h(T-r) dr 


If the use of an integrate and dump filter is assumed, ie. hit}=1 for te (-1/2,7/2), then 


(3.22) reduces to 
A ( 7) hy neq2 2 75,2 


ae 2 15,2 
te t7tt) 20% 4 Pret 20) ar (3.24) 
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Also noting that 
Erf(u) = (21) fo" exp(-x®) dx (3.25) 
Equation (3-24) reduces to: 


An =(m/hbv) (qRA)” <G?>( 27) */o. 
[Pyerf(T/(2véc)) + P,ferf(3T/(2véc)) — erf(T/(2rb0)){] (326) 


In a similar manner one can derive the n'th cumulant for the case of intersymbol 
interference, when the highest level signal occurs between time slots occupied by the 


lowest level: 
dn=(n/hv) (qRA)® <G®>( 27) '/2o. 
[Pyert(T/(2véc)) + Piferf(3T/(2véc)) — ert(T/(2v60))}] (327) 


Equations (3.14)-(3.20), (3.26) and (3.27) can be solved to find the optimal 
thresholds and signal levels for a fixed error probability. Conversely, given the signal 


levels and thresholds, one can calculate the error probability. 


3.0.2 Numerical Example 


We proceed here, with the calculations of signal levels and thresholds needed to 
achieve an error probability lying in the range of 10-* to 10-}° The system parameters 
used in this example are given in Table 3.1. The output optical power waveshape, git), is 
assumed to be Gaussian with variance,t,. 

Equations (3.14)-(3.20) and (2.86) can be used to calculate the optimum signal 
levels and thresholds for different error probabilities. The results obtained from solving 
these equations are given in Table 3.2. These were calculated assuming no intersymbol 
interference (i.e. the effect of git) outside the time slot under consideration is ignored). 
From Table 3.2 one can see that the optimum signal levels and thresholds are not equally 
spaced. This results from the fact that the shot noise term increases with increased signal 
level, thereby making it necessary to keep the thresholds closer to the lower levels (so as 


to maintain areas A,, A,j,.....A,, shown in Fig 3.2, equal). 
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TABLE 3.1 


PARAMETERS USED IN NUMERICAL CALCULATIONS 


DEFINITION 


Quantum efficiency of the APD 


Amplifier inpur resistance 


Amplifier gain 


Average avalanche gain 


Optical carrier frequency 


APD Ionisation ratio 


Variance of Gaussian optical pulse edge 


Time slot available per pulse 


Thermal noise 


VALUE 


700 ohms. 


1500 


10 


3.614x10" Hz. 
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The results obtained above change drastically when one takes into account the 
effect of intersymbol interference. Fig. 3.3 shows graphically the worst case intersymbol 
interference for a 4—level system, assuming Gaussian shaped optical pulses at the input 
of the detector. The worst case intersymbol interference occurs when the lowest level 
signal is flanked on adjacent time slots by the highest level signals. Fig.3.3 clearly shows 
that the effect of the worst case intersymbol interference is to increase the amplitude of 
the lowest level signal. There are therefore greater chances of the level—1 signal being 
detected as the level-2 signal. As a consequence of this, the overall error probability 
increases. The increase in error probability due to intersymbol! interference can also be 
shown analyticaliy. Equations (3.26) and (3.27) give the cumulants of the output voltage 
assuming intersymboi interference. Intersymbol interference is seen to increase the 
cumulants of the output voltage. As a result of the increase in the mean value and 
variance of the output voltage, the areas A, and A,, representing the conditional error 
probabilities P[ 1/2} and P[2/1], increase. The overall error probability increases due to the 
increase in the conditional error probabilities. 

Using (3.27) along with (3.14) to (3.20) and assuming the same system parameters 
as before, one can recalculate the error probability taking into account the intersymbol 
interference. It is seen that, for an optical pulse with a variance, t, , equa! to one-fifth of 
the time slot, the error probability increases from 10°* to €75 x 10° The error rate 
increases rapidly for larger values of t,.. For an optical pulse with t equal to two-fifths 
of the time-slot the error probability deteriorates to a value of 0.25. Such a condition is, 
however, an extreme case of intersymbol interference. Appendix 3 gives programs to 
calculate error probability for various values of t.. 

A plot of error probability as a function of the mean avalanche gain, <G>, of the 
APD is given in Fig. 3.4. The value of <G>is varied from 10 to 200. It is seen that the 
error probability decreases as <G> increases up to a vaiue near 100. Above <G> = 100, 
the error probability starts to increase again. This behaviour is due to the fact that for 
lower values of <G>, the system is themal noise limited However, for <G> greater than 
100 the system starts to become shot noise limited. The optimum gain for the above 
system is therefore near 100. This sort of behaviour has also been observed for 


two-level systems [4]. 
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(a) The effect of the worst-case intersymbol interference (ISI) on the pulse 


shape. The dotted line represents the pulse shape without ISI and the 


thick line represents the pulse shape after ISI. Note the increase in 


the mean amplitude of the lowest level. 


x* 


(b) The effect of the worst-case ISI on the PDF of the output voltage. The 


mean level and variance of ‘signal-1' is larger than that without ISI. 


Figure 3.3 The effect of the worst-case intersymbol interference on the pulse shape and 
the probability density function due to the highest level signal flanking the 


lowest leve! signal on both sides. 
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Figure 3.4 Plot of error probability as a function of the mean avalanche gain for a 4-level 


PAM system. 
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Fig. 3.5 shows the variation of the error probability as a function of the average 
optical signal, Pavg at the input to the detector. The error probability is seen to decrease 
with increased optical power reaching the detector. These results were observed for the 
2-level case by Mansuripur et. a/. [19]. Fig. 3.5 also shows the variation of error 
probability as a function of the average optical power reaching the detector, taking into 
account intersymbol interference. The variance of the pulse, t_ is assumed to be 
one-fifth of the time slot. The programs used to generate these plots are given in 


Appendix 3. 


3.1 EXPERIMENTAL SETUP FOR THE 4-LEVEL PAM SYSTEM 

A complete four-level PAM system was designed, built and tested in the 
laboratory. A major objective for the experimental part of the project was to evaluate 
the actual performance of the 4-level system and to determine what practical advantages 
this approach has over other systems. A further objective of this experiment was to was 
to verify the validity of the Gram—Charlier series for the characterisation of the PDF of 
multilevel systems. 

With the above objectives in mind, the following measurements were made using 
the experimental 4-level PAM system: 
a Error probability as a function of the average optical power at the receiver input. 
2. ‘Error probability as a function of the average avalanche gain of the APD. 
The remainder of this chapter describes in detail the circuits used for the experiment, as 


well as the experimental results. 


3.1.1 Overall Experimental Setup 
In Fig. 3.6 the 4-level PAM system and test equipment that were used for the 


experiments are shown. The 1645A Data Error Analyser (manufactured by Hewlett 
Packard) was used to generate a pseudo-random binary sequence (PRBS) at a rate 
determined by an external clock. For this experiment, the external clock was set at 4 
Mb/sec, which is close to the maximum usable bit rate of the 1645A Data Error Analyser. 


The fiber-optic system is, however, capable of operating at much higher rates. 
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Figure 3.5 Plot of error probability as a function of the mean optical power reaching the 


detector. 
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Figure 3.6 The experimental four-level PAM fiber optic system. 
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The PRBS is applied to a 4-level modulator which in turn intensity modulates the 
LED output by means of a driver circuit. The optical signal is coupled to a fiber of 600m. 
length. The output of the fiber is coupled to an APD which coriverts the optical signal 
back into an electrical signal. The received signal is then routed through a pre-amplifier 
and main amplifier. The amplified signals are fed to an integrate-and-dump filter to 
eliminate some of the noise. Finally the signal is demodulated and fed back to the Data 
Error Analyser as the data input. The Data Error Analyser compares the system input 
sequence with the system output sequence and displays the BER. 


Each of these stages are described in detail in the following pages. 


3.1.2 The Data Error Analyser 


The HP 1645A Data Error Analyser is well suited for data communication systems 
and has the capability of simultaneously measuring the following parameters: bit error 
rate (BER), block error rate (BKER), clock slip rate, carrier loss rate, skew, jitter and total 
peak distortion. These parameters are measured during a test sequence, digitally stored, 
and read out later. ) 

For the purposes of this project, the BER and BKER measurements were of 
greatest interest, although the other measurements did provide added insight into the 
operation of the system. 

The 1645A computes the error rate by correlating the received bit pattern with a 
local replica of the transmitted data sequence. A pseudo-random binary test sequence 
(PRBS) can be generated of length 2?°-1, 24-1, 2°-1 or 2-1 bits. Also sequences of 
type 7:1 (7 zeroes followed by a one), 3:1 and 1:1 can be generated. The 1645A divides 
the test data sequence into 1000-bit blocks for test analysis. The BER is the number of 
times a bit error occurs during the full test sequence. The BKER is an indication of the 
number of received blocks containing one or more errors. The BER test is used to 
measure the overall link quality. BKER when compared to BER on the same data base, 
provides an indication of the burst characteristics of the errors. 

If errors occur in a distributed fashion (high BKER) then other tests such as skew, 
jitter or total peak distortion should be reviewed. Also high BER with low BKER indicates 


that phase hits (ie. clock slips) and dropouts should be checked. The SKEW measurement 
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indicate the percentage of ones changed to zeros. Ideal skew characteristics of 
erroneous data is 50% ie. errors should be equally distributed between ones being 
mistaken for zeros and vice-versa. If errors are skewed, there are two possible causes; 
errors are being repeated (patterned) or a decision threshold is set at an incorrect level, 
so that errors are weighted in one direction or the other. 

The 1645A also has a digitally matched filter which can be switched in or out 
ahead of the error detection circuit. The filter removes the high-frequency noise and 
prevents data sampling on spurious pulses. 

The flowchart of Fig. 3.7 summarises the measurement troubleshooting 


procedure [32]. 


3.1.3 Four-level Modulator and LED Driver 

The 2-level PRBS signal coming out of the 1645A Data Error Analyser is 
converted to a 4-level signal using the circuit shown in Fig. 3.8. The basic function of this 
circuit is to combine two binary pulses and generate one out of 4 possible signal levels 
based on this sequence. The possible sequences are; 00, 01, 10 and 11. A sequence of 
OO produces the lowest of the four signal levels and a sequence of 11 produces the 
lar gest. 

U;, and U,, form a 2-bit shift register. After 2 clock pulses the inputs to U,., 
and U,, contain the first and second bits respectively. The outputs of U,, and U,, are 
coupled to transistors Q, and Q, on the application of a clock pulse from U, via U,. Since 
Q, and Q, are emitter coupled, the current flowing through the LED is the sum of the 
individual emitter currents which depend on the individual base voltages which in turn 
depends on the signal sequence. The choice of the optical source was based on 
considerations of reliability, linearity, dynamic range and cost. The Northern Telecom 
NT-40-3-30-3 LED was found to best fit these requirements and was consequently 


selected as the optical source. 
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Figure 3.7 Algorithm for the troubleshooting of the optical fiber system. 
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Figure 3.8 Circuit diagram of the 4—-level modulator and LED driver. 
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3.1.4 The Choice of the Fiber 


A typical fiber-optic system design requires the selection of the link length, the 
repeater spacing, the rise time budget and the loss budget. The aim of the present 
experiment was to observe the effect of the various noise sources on the BER. A major 
consideration in choosing the fiber type was that it should have a sufficiently narrow 
bandwidth that intersymbol interference will be observed. To observe intersymbol 
interference for a binary sequence transmitted at 4 Mb/sec., the bandwidth of the fiber 
should be lower than 8 Mhz. [33]. A GI Fiber (Siecor Super Fat Fiber Cable 155) with 
bandwidth of 5 Mhz.km. was chosen. This fiber has a very large attenuation of 35 db/km. 
Keeping in mind the sensitivity of the receiver, the length of the fiber was chosen to be 
600m. (these fibers are drawn to lengths of 100m., 400m. and 600m). Other 
characteristics of this Siecor fiber are its comparatively high numerical aperture (i.e. its 


light gathering ability) of 0.4, its high strength and its flexibility. 


3.1.5 Receiver Section 

The front-end of an optical receiver consists of a photodetector (generally a PIN 
diode or an APD) along with some form of amplifying stage or stages. The overall 
combination of the detector, amplifier and subsequent filtering is designed to respond to 
the input light signal in such a way as to provide an output pulse shape (usually with raised 
cosine spectrum) appropriate for presentation to the digital decision circuit. The 
combination, referred to as a linear channel, must therefore have sufficient bandwidth to 
respond properly to the input pulse. It should also contribute as little noise as possible in 
order to give good optical sensitivity. 

There are two approaches normally used for receiver design, i.e. the high input 
impedance method and the feedback method. In the high input impedance design, note is 
made of the fact that bit error probability is reduced by increasing the input resistance 
and decreasing the input capacitance. This, however, results in the integration of the 
detected signal. To compensate for this an equaliser stage is necessary at the end of the 
amplifier stage. Quite often the equaliser is just a differentiator. The equaliser does not 
introduce much noise because it follows the amplification stages. Most of the noise 


contribution occurs at the input. Thus a high-impedance equalised amplifier is capable of 
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extremely low noise while retaining the pulse shape. There are, however, a few problems 
with the high-impedance approach. The transfer function of the equaliser depends on the 
value of Ry and Cy (see Fig 24,R_ =R, // Ry andC, =Ca + Cy ), which depend 
in turn on parasitic capacitances and on the 6 of the transistor. As a result, individual 
amplifiers may have to be individually equalised. In addition, temperature compensation 
schemes may be neccessary. This complicates the equaliser design. Another problem is 
the reduction of the dynamic range due to the build-up of charge on the amplifier input 
Capacitance. This problem can be eliminated by actively discharging the input capacitance, 
but this introduces more complexity into the receiver design. 

Another preamplifier type, commonly employed for front end design is the 
transimpedance or current-to-voltage converter. The transimpedance amplifier has wider 
bandwidth (BW), greater dynamic range and a noise performance almost as good as the 
high-impedance amplifier [34]. The increased bandwidth is due to the feedback scheme 
being used. The wide dynamic range arises because the low-frequency signal has a 
closed loop amplification (which is controlled by external elements). It is the low 
frequency signal which is increased by front-end integration thereby causing the 
amplifier to saturate prior to the equaliser. Improvement in the dynamic range is 
proportional to the ratio of the open-loop gain to the closed loop gain. The major 
drawback of the transimpedance approach is the fact that it tends to oscillate for high 
gains. 

The feature of the transimpedance amplifier that makes it desirable for the 
present use is that, compared to an unequalised amplifier not employing feedback, it is 
less noisy for a given BW (the resistor noise is reduced) or alternatively has more BW 
for a given noise level. It does this at the expense of extra optical power (typically 1 
db.) to achieve a given error rate. Circuit simplicity, eliminating the need to employ 
equalisation and obtaining increased dynamic range were judged worth the 1 db. loss in 
sensitivity. 

Fig. 3.9 shows the circuit diagram of the FET preamplifier that was used for this 
experiment. The circuit used here is a modification of a GaAs FET amplifier designed by 
Ogawa and Chinook [35]. The design of the first stage requires particular attention since 


it normally contributes most of the amplifier noise. Since this particular feedback 
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Figure 3.9 Circuit diagram of the FET preamplifier. 
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configuration is considered for application below 100 Mhz., a junction FET was chosen 
as the active device. At this frequency, commercially available FET’s are generally less 
noisy than commercially available bipolar transistors. 

The first stage FET is operated optimally in the common source mode. In this 
mode, the input impedance is high, as dictated by (2.31). A common-base stage follows 
the F.E.T. Such a configuration used with the a bipolar transistor minimises the Miller 
Capacitances and produces a high open loop gain. 

The emitter follower stage, using an n-p-n transistor, is used after the 
common-base stage. The feedback resistance is connected from the emitter of this 
transistor to the base of the FET. The value of the feedback resistor was chosen to be 
equal to 20 k_ in order to avoid oscillation of the preamplifier. Further amplification is 
provided by the next stage. The final stage is an emitter follower which is used because 
of its low output impedance. A low output impedance is necessary for impedance 
matching with the main amplifier (Keithley 105 Pulse Amplifier) that was connected after 
the preamplifier. The voltage gain of the preamplifier was measured to be 30 and the 


gain of the main amplifier was set to 100. 


3.1.6 Integrate-and-Dump Filter 

The integrate and dump filter, at the output of the receiver main amplifier, helps 
to smooth out the noisy amplified voltage and thereby increase the signal to noise ratio. It 
is also needed in order to prevent oscillation of the comparator circuits used in the 
demodulator. 

Fig. 3.10 is a circuit diagram of the integrate and dump filter that was used for the 
4-level PAM communication system. The basic requirements of this filter are a slow 
integration time (between 100 to 200 ns.) and a very quick dump time (less than a few 
nanoseconds). To do this, an RC integrator has been used with a diode connected in 
series with R to avoid discharge of the capacitor. A VMOSFET switch (VN10KM) is 
connected across the capacitor. Whenever the gate signal is above the pinch-off voltage 
of the FET (in this case 2.0 V), the FET is turned on such that the drain-source resistance 


(r ) is very small (approx. 5 ohms). 
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Figure 3.10 Circuit diagram of the integrate-and-dump filter used for the 4—level PAM 
system. Also shown are the input and output waveforms and the control 
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3.1.7 Demodulator 


The output of the integrate and dump filter is a sequence of 4-level signals which 
has to be converted back into a 2-level sequence before any error rate measurements 
can be made. This conversion is achieved by the demodulator circuit. The circuit diagram 
of the demodulator that was designed and built for this project, is given in Fig. 3.11. The 
clock signal to the demodulator is the same 4 Mb/sec. clock that was fed to the 
modulator, except that it is delayed by a certain time period. The delay of the clock 
should be set so as to match the delay in the transmission of the signal through the fiber. 
Since this delay is not known exactly it was adjusted so as to give the lowest value of 
BER. 

The first stage consists essentially of three comparator circuits (LM361) which 
decide which of the 4 levels was transmitted. The rest of the circuit consists of NAND 
and OR gates which reconstruct the 2 level signal based on the output of the 


comparators. Fig 3.12 presents the timing diagram for this circuit. 


3.1.8 Error Rate Measurements 

The demodulator output is fed to the data input terminal of the Data Error 
Analyser. The Data Error Analyser correlates the received bit pattern with the transmitted 
data sequence and determines the number of errors occuring in a certain amount of data. 
It automatically computes the BER for that length of data. 

Measurements were made of the BER as a function of the average avalanche gain 
of the APD and as a function of the average optical power reaching the detector. The 
results obtained have been plotted in Fig. 3.13 and in Fig. 3.14. The measured signal levels 


and error probabilities have been given in Table 3.3. 


3.2 DISCUSSION OF RESULTS 
In order to compare the experimental results with theoretical calculations it is 
necessary to determine the following system parameters: 
1. | Thermal noise of the amplifier. 
2. The value of the ionisation coefficient, K. 


3. The value of the mean avalanche gain as a function of the reverse bias voltage. 
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Figure 3.11 The 4-level demodulator circuit diagram. 
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Figure 3.12 Timing diagram of the 4-level PAM demodulator. 
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Figure 3.13 Error probability as a function of the mean avalanche gain of the APD for an 


experimental 4 level PAM system. 
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Figure 3.14 Error probability as a function of the average optical power at the input to 


the receiver for a 4 level PAM system. 
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TABLE 3.3 


MEASURED VALUES OF SIGNAL LEVELS AND THRESHOLDS FOR THE 
EXPERIMENTAL FOUR-LEVEL PAM SYSTEM 
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4. The quantum efficiency of the APD. 

The thermal noise was determined by measuring the output signal when there was no 
input signal to the amplifier. The value of thermal noise was estimated to be 15 mV, from 
measurements made with the Spectrum Analyser. The value of the ionisation coefficient 
is highly dependent on the individual device. For Silicon, its value is typically in the range 
0.01 to 0.04. The calculated error probability did not vary much over this range of K, and 
a value of 0.02 was chosen in all the theoretical calculations’ 

The value of the average avalanche gain as a function of the reverse bias voltage 
for an APD is determined metetrally by fabricating a PIN diode and an APD on the same 
chip and taking the ratio of the output currents. This ratio directly gives the average 
avalanche gain of the APD. A simpler method, but not as accurate, is to determine the 
reverse bias voltage beyond which one actually observes increased output signal. This 
reverse bias voltage, or any smaller voltage, then can be assumed to correspond to a 
unity avalanche gain. For any other setting of reverse bias, the average avalanche gain can 
be assumed to be the ratio of the output signal for that setting to the output signal for an 
average gain equal to one. This method is accurate to within 10% °. Fig. 3.15 is a plot of 
the average avalanche gain of the APD as a function of the reverse-bias voltage. 

The value of the quantum efficiency of the APD was not supplied by the 
manufacturer. The value of the quantum efficiency can vary typically from 0.6 to 0.9. The 
error rate was observed to be fairly insensitive to the value of the quantum efficiency 
over this range and was arbitarily chosen as 0.85 for all the theoretical calculations. 

Figs. 3.13 and 3.14 show the variation of the error probability as a function of 
the average avalanche gain of the APD and as a function of the average optical power 
reaching the detector. The experimental results are Compared to the results obtained 
theoretically for a range of values of the thermal noise, th: It is seen in this graph that 
the error probability is highly sensitive to the value chosen for Cth: 

The effect of intersymbol interference on error probability is also shown in Figs. 
3.13 and 3.14 by considering different values of rise (or fall) times, t, , of the optical 


pulse reaching the detector. In actual practice, it is difficult to determine t, accurately. 


3Private Communication, R.J. Mcintyre, RCA labs, Montreal. 
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Figure 3.15 Mean avalanche gain of an APD as a function of the reverse bias voltage 
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Other factors which could affect the results are the inconstant delay of the clock 
pulse to the demodulator. The theoretical results were obtained assuming that the 
decision is made at the same position during each time slot It was observed, however, 
that the delay of the clock pulse to the demodulator tended to vary with time. In addition, 
the readings of the optical power were accurate to ++1 dB. It has already been shown that 
the error probability is highly sensitive to changes in the received optical power. From 
Fig. 3.5, even a 3 db. change in power can cause the error probability to change by a few 
orders of magnitude. Thus the errors in reading the optical power could result in a 
mismatch between experimental and theoretical results. 

Despite the various inaccuracies that could result from all the above factors, the 
results shown in Figs. 3.13 and 3.14 match fairly well with the values obtained using the 
Gram-Charlier series. As seen in Fig. 3.13, there is a slight discrepancy between the 
theoretical and experimental values of error probability at higher values of avalanche gain 
(i.€., above 10). The discrepancy arises because the effect of saturation of the main 
amplifier, which occurs at high avalanche gains, was neglected in the theoretical 


calculations. 
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CHAPTER IV 
MULTILEVEL PULSE WIDTH MODULATION 


The advantages of multilevel signal transmission over regular binary signal 
transmission has been discussed in Chapter Ill. The analysis thus far has, however, been 
restricted to PAM systems. Multilevel PWM transmission is shown, in this chapter, to be 
an attractive alternative to multilevel PAM. 

PWM and PPM transmission have been recognised to be more immune to noise 
than PAM transmission [37], provided the transmission channel has a sufficiently large 
bandwidth. PAM is likely to have fluctuations in the output amplitude due to the various 
noise sources in the system. Such fluctuations in the signal amplitude do not, however, 
affect the signal content of PWM or PPM systems. This is due to the fact that PWM and 
PPM techniques involve only the transmission of timing information and not amplitude 
information. PWM and PPM can, therefore, fully utilise the wideband characteristics of 
optical fibers in order to decrease the bit error rate. 

In this chapter, the pulse widths and threshold levels required for a near—optimum 
4-level PWM system are determined. This extends the work of Personick [12] for binary 
systems and Hullet and Muoi [31] for multilevel PAM systems. It is shown that, given a 
PDF for the noise amplitude, the PDF for pulse-width can be computed. Thus, the 
analysis of error probability for PWM can proceed in the same manner as that for PAM. 

Two different PWM schemes, one employing a low pass filter (LPF) and the other 
an integrate and dump filter (IDF), are analysed. It is shown that the LPF scheme can 
achieve very low error probabilities, provided that a wideband fiber is used. The IDF 
scheme is shown to be advantageous in systems where such wideband fibers are not 
available. Examples comparing the two different schemes are given in this chapter. This 
analysis has been restricted to Gaussian shaped pulse edges at the fiber output. Also, for 
the sake of simplicity, a Gaussian PDF for the overall noise has been assumed. This 


analysis could be extended to a non-Gaussian case using the Gram—Charlier series [19]. 
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4.1 MATHEMATICAL MODEL OF THE MULTILEVEL PWM SYSTEM 


The fiber optical system model is shown schematically in Fig. 4.1. The message 


source generates a sequence of Statistically independent random variables a; with 


J 
values {1,2,3,4} : 


oo 
383) ) a, SCtaiT) 4 ee Cle 2384) (4.1) 
i= 
where T is the time slot available for each pulse. 


oo 


The encoder converts: the source signal into a symbol! sequence suitable for transmission 


as PWM: 


foe) 


x(t) =] r (s,, t-iT) 4 cl (me2g3s 4) 


1 = —& 


(4.2) 


where r(s,t) is an ideal rectangular pulse of duration s, 
and we assume that: s, < s, <s, < s, < T. 

The encoded pulses, x ,{(t), suitably amplified, drive the optical source. The optical 
pulses generated by the source will depend on the bias current, the modulating signal 


Current and on the source response. We denote the output of the optical source as: 


ay (t) 2 ) Soy ee t-iT) jt cs 23 Sh, 4) (4.3) 
ah = ==O0 


where r,(t) represents the pulse shape at the output of the optical source. 
The optical signal suffers attenuation and distortion along the fiber channel. 
Assuming that the fiber has an impulse response N¢ (t), then the optical power, Pit), at the 


output of the fiber is: 


P (t) = he (t) * Ty (s;> t) Sec Clem ord) (4.4) 


The fiber response, h¢ (t), is assumed to be such that P(t) has Gaussian leading and falling 
edges (as shown in Fig. 4.2). This assumption has been shown to be fairly accurate, under 
certain conditions [38]. For purposes of analysis, the pulse is assumed to start and end at 
the point where the amplitude of the Gaussian rising and falling edges is 10% of the 
maximum pulse height. Neglecting the tails of the Gaussian pulses (i.e, the areas below the 
10% level) does not affect the calculated error rates for the LPF method of detection and 


causes only a slight error in the results for the IDF scheme. Thus the optical pulse shape 
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Figure 4.1 Block diagram of a multilevel PWM optical-fiber communication system. 


Figure 4.2 Assumed pulse-shape of the optical output powerof the fiber. 
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at the output of the filter, within the i’th time interval, can be represented as: 


mata iTe= Ee 
a ateteemoiea me 
2t 2 


r exp 


_ 


for tc (iT, iT+2.144t_) 


D (Gee shay & P. for tlew (it +2.144t WI Sriee +0.144t J 


m 


2 
' = Sf -_ — ' 
P om (t - iT t. 0.144t' ) 
2t1 2 


for tc (iT +t, +0.144t , iT +t, +2.288t,) 
(4.5) 


where: io is the time variance of the Gaussian edges and 2.144t0 is the rise time of the 
opticai pulses (from 1/10 th of the max. value to the max. value). 

The receiver section consists of an avalanche photodiode (APD), an amplifier, a 
filter and a decision circuit (Fig. 4.1). The output current of the photodetector can be 


described as a filtered randomly multiplied Poisson process. The current from the APD is: 


k(tp) 
i(t) = } &, dh, (t - 
eae ee *n? (4.6) 


where q is the electronic charge, h, (t) is the response of the APD and where the 
number of primary electrons k(t ,), the random time points t , and the random avalanche 
gain g,, are statistically independent random variables. The amplifier section introduces 
Gaussian thermal noise. 


The average voltage at the filter output (Fig. 4.1) at time t is given by (2.78) 
E 


vec) a= oS (qRA) a) Ie Gal (Qe G3) ce nad 


where hit) is the ideal impulse response of the filter (Fig 4.3a & 4.3b) 

The decision circuit for the LPF scheme can consist of a threshold detector 
(slicer), an integrator and a sampler. The output of the slicer is a “clean” pulse with sharp 
leading and trailing edges. The slicer preserves the width of the received pulses (as 


modified by noise). In other words, the slicer eliminates all the amplitude noise except 
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Figure 4.3 Impulse response of (a) an ideal IDF and (b) an ideal LPF. 
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that in the neighbourhood of the threshold, thus greatly reducing the final output noise of 
the system. One method of measuring the pulse width is to integrate and sample the 
sliced output every T seconds, yielding a PAM wave. The IDF also produces PAM signals. 


The original pulse samples can then be regenerated for both the LPF and the IDF 


schemes. 


4.2 NOISE IN A 4-LEVEL PWM SYSTEM 
The analysis of the noise properties of a 4-level PWM system are considered 
here. This analysis has been performed for both the LPF and the IDF schemes. These 


results are easily extendable to any multilevel PW/M system. 


4.2.1 LPF system: 
Fig. 4.4 shows the relationship between the noise amplitude and the timing error 
at the output of the LPF. The relationship between AT (i.e. the timing error due to noise 


amplitude n, at time t) andn4 is: 
ob ay a 
AT ale (4.8) 


The threshold level Ath is chosen so as to minimise the effect of n, on AT i.e., 


ve dv (t) yt Veen 
ce hry; ee — (4.9) 
dt max 
Ae t ve 
fo) 
where: 

t,- is the time variance of the Gaussian—shaped pulse edges at the output of the 
ete 


Vmax is the maximum voltage at the output of the LPF given by : 


t. + 0.144t) + aT, 


i 
= (pels ' n 
Lae ( hy) (qRA) <G>f P(t) h(t, + O.144t' + iT - t)dt 
(4.10) 
ath 
The normalised random output voltage, at time t, with mean m and variance q is: 
t] 
X -m Me 


a ve (4.1 1) 


tarwitiatoa ste mareye MWS (ovel-O's Fora 
oanaT nanesiga SO) edn Bria V5) sett a pt 


ore gnimit st brie Souilicena seve sit hseiltiad Giendipivy ent aworte hh 
salon ah suhh ioe HA Sat au he rotted cleanablitee nA? 7.) att lo tus 
a gi 'Gra omit ta yr at 


ie [ 1 = ie i. 7 a. 


Ai GE ite fs98 to aT eli ata core Poll aves © v 7 
ry 2 = 
y ee a 


ae! . aie 
ait Yo Sleue: oed 26 aggae “elug sequne rapa WE To Sorkeaw amit wit a ee) Sa 


qa neVig V4) ect ta tec Sey) : a un eet 
TY 43 


; : Beta a ; ae 
apg Ter Seid. + old Ga Yam tate) GE) © 2a 


crab 


[ . 4 - 
a ada ial eli 1 ait tw ig 
| - 


$3 


Magnified View 


Slicing Level 
Slicing Level 


— — — — Signal Without Noise 


Signal Plus Noise 


Figure 4.4 Relationship between the noise amplitude and the timing error at the output of 


the LPF. 
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or 


t t (4.12) 


where: 
X' is the non-normalised random voltage and 
X is the normalised random voltage. 


The variance, oj , is given by (3.20): 


i 
we SL aS /2 


where Ae is the second cumulant of the random voltage at the output of the filter and 
Oy, 1s the thermal noise of the amplifier. 


ty + 1.144t" + iT 


N= GL) Gaye <C7> é ' 
2 hy q PCr) eh (t, + 1.144e' +47 - 1) dt (444) 
ES 
<G’> is the second moment of the gain of the APD given by [39]: 
<G2> = Shak 2 3 
al -k)@ <C> + 2(1-k) <G>se4a <cs (4.15) 
where k is the ionisation coefficient of the APD. 
Substituting (4.8) into (4.11) we get 
t veo 
t X 
At =————— (4.16) 
Vv 
max 
or 
Le eel (4.17) 
where K is a constant for a particular pulse-width, given by 
yaa 
KES (4.18) 
Omi tiny £6! 
ET Yo} 


From (4.16) the final relationship between the PDF of the normalised noise 
amplitude X and the PDF of the pulse width At , is derived as: 


P(At) = P(X (At)) eee (4.19) 
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Assuming that the PDF of X (the normalised random variable) is Gaussian, we have: 


i 
— -X 
09) Ie E30 eee, (4.20) 


If t represents the random pulse width, then; 


ito a PY (4.21) 


From (4.17) and (4.19) the PDF of Tt is: 
K 


XE ayo ™ oe Pe {kK Cue aur 


2 


(4.22) 


where t; 1s the average pulse width. Fig. 4.5 shows the PDF for a 4-level PWM system 

using the LPF scheme. Here: 
t, (i=1,2,3,4) correspond to the mean pulse widths. 
Tj (j=1,2,3) correspond to the optimal thresholds. 


The above theory can be used to compute the error probability as shown in Section 4.3. 


4.2.2 IDF system: 


An alternative to the LPF system of detecting pulse width is the IDF scheme. The 
pulses at the output of the preamplifier are integrated and a decision is made about the 
pulse width, depending on the amplitude of the integrated output. 


The PDF of the random output voltage is assumed to be Gaussian, i.e. 


60 eae 
1 : 
fe eee FP eek 
2T O. Oo. (4.23) 
i i 


where ora the noise variance term for the ith level. 

The procedure for calculating GIs the same as for the LPF scheme, except that 
in (4.12) and (4.13) one must use the impulse response of the integrating filter, shown in 
Fig. 4.3b. 


Fig. 4.6 shows the PDFs for a four-level PWM system using the IDF scheme. 
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Figure 4.5 Probability density function of the random timing error at the output of a 


four-level PVWVM system using an LPF. 


a7 


Figure 4.6 Probability density function of the random timing error at the output of a 


four-level PWM system using an IDF. 
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xX" is the random voltage at the output of the IDF. 

m; is the mean voltage at the output of the IDF corresponding to a pulse of width 
_t,; being transmitted. (i= 1,2,3,4) 

ue is the amplitude threshold (chosen using the equal area criteria). These 


correspond to the time thresholds GF 1.253) 


The error probability can then be computed as shown in the next section. 


4.3 ERROR PROBABILITY FORMULATION FOR A 4-LEVEL PWM SYSTEM 


The probability of error P(€ ), for the 4-level system is: 


P(e) = Py P21 ee P» (P(1/2) + P(3]2)) P (P(2|3) + 


+ P(4|3)) + P, (P(3]4) (4.24) 

where: 
P(ilj) = probability of receiving level i when level j was transmitted. (i,j = 1,2,3,4) 
andp,; (i= 1,2,3,4) is the probability of sending a pulse of widtht; . Assuming that the 


signals have been efficiently randomised so that each of the four levels is equally 
probable, then p will be equal to 1/4. Secondary error terms such as P(3!1) and P(411) 
have been neglected in (4.23). 

From Figs. 4.5 and 4.6 it is noted that the above conditional probabilities 
correspond to the areas A,B,C,D,E and F. Using the equal area criteria and (4.23), each of 
these areas must be equal to two-thirds of the overall error probability. The error 
probabilities can be found from (4.21) and (4.22) for the LPF and IDF schemes, 
respectively. 


Thus; for a LPF: 


co K) 2 
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si 


The other areas B,C,D,E and F can be derived in a similar fashion. The following relation is 


used to compute the optimal pulse width levels t, once the m,’s have been calculated. 
t, + 2.288 t' 
OG 


m, = Gree) (qRA) <G> P(t)h(t, ae PAPAS Lae E)edt 


ie) 


(4.26) 


where hit) is the impulse response of the IDF as shown in Fig. 4.3b. 
Thus, given the pulse levels and thresholds, the error probability can be computed 
using the above equations. Conversely, given the required error probability, the 


near-optimal pulse width levels and thresholds can also be computed. 


4.4 NUMERICAL EXAMPLES 

We proceed here with the calculations of pulse-width levels and thresholds 
needed to achieve an error probability lying in the range 10-* to 10°° The system 
parameters chosen for this example are given in Table 4.1. 

Unlike PAM, where pulses of unlimited amplitude can theoretically be used, in 
PWM the maximum pulse width is equal to the time slot. The pulse amplitude can, 
however, be set to achieve any particular error rate. For this example, the maximum pulse 
amplitude is assumed to be constant at 15 nW, with a minimum pulse width, t,, set to be 
10 ns. Using equations (4.24a) to (4.24b) (and other equations derived similarly) one may 
obtain the thresholds and levels for the LPF and the IDF schemes for error probabilities 
varying from 10-4 to 10°!*. These are given in Table 4.2. 

A graph of error probability vs. P,, 1s shown in Fig. 4.7 for both the LPF and IDF 
schemes. A mean avalanche gain of 100 has been used. The results for 3 different rise 
times have been shown. Fig. 4.8 is a plot of error probability vs. <G> for the LPF and IDF 


schemes. These plots have been made for different values of PR, (SnW, 15nW and 
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TABLE 4.1 


PARAMETERS USED IN NUMERICAL CALCULATIONS 


: 
. = 


<G> Average avalanche gain 100 


4 
v Optical carrier frequency 3.614x10 Hz. 


K APD Ionisation ratio 
ty Variance of Gaussian optical pulse edge 


Sth Thermal noise 3 mV 


Ww Bandwidth of the LPF ™x10° rad/s. 


TABLE 4.2 


REQUIRED PULSE WIDTH LEVELS, USING AN IDEAL LPF AND _ IDF 


FOR A RANGE OF ERROR PROBABILITIES (TIME IN nanoseconds 
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Figure 4.7 Error probability versus F,, for both the LPF and IDF schemes. 
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25nW). Finally, for the LPF scheme, error probability is plotted as a function of the rise 
time (t, ) of the optical pulse edges at the input of the receiver in Fig. 4.9, for different 
power levels. Figs. 7,8 and 9 have been plotted using the same values of pulse widths as 


have been calculated for the LPF scheme above. 


4.5 DISCUSSION OF RESULTS 

From Table 4.2 it is seen that the error rate decreases rapidly with increased slot 
time. A relatively small change in t, , from 28 ns. to 31 ns., reduces the P(é ) from 10-7 to 
10°°. From this table it appears that the LPF scheme gives lower P(€ ), for a given slot 
time, than the IDF scheme. However, this is only true for the particular rise time chosen 
(t_= 1 ns). For higher values of rise time, the IDF scheme becomes better than the LPF 
scheme. 

This effect is seen clearly in Fig. 4.7, where, although the value of t, has very 
little effect on error rates in the IDF scheme, it is a crucial factor in the LPF method. 
Using an input power P,, = 15 nW., the P( € ) for fs IenSaiSe lO evwnereasstOmte is 
ns. the P(é)"is—l0-'", 

From Fig. 4.8 it is observed that, for P,, = 15 nW., the error probability for the 
IDF scheme initially decreases as gain is increased. For gains higher than about 90, the 
P(€) starts to increase as gain is increased. Thus, there exists an optimal gain for the IDF 
scheme. For different values of P,,, , one finds a similar behaviour except that the optimal 
gain decreases with increasing P,,. 

A similar effect is observed theoretically for the LPF scheme, except that here 
the optimal gain is near 10. If the thermal noise term is increased, the optimal gain of the 
LPF shifts to higher values of <G>. 

From Fig. 4.9, which is a plot of error probability vs. t, for a LPF scheme, it is 
clearly seen that P( € ) decreases greatly as t , decreases, ie. as the system bandwidth 
improves. Theoretically, one can reduce P(E ) to any desired level provided that the 
bandwidth of the system is large enough. The LPF scheme has two degrees of freedom 
in reducing P(E) viz. either by increasing power levels or by increasing the system 
bandwidth. This is not possible for the integrate and dump scheme where the only way to 


reduce P(€) is by increasing the power level. 
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Figure 4.9 Error probability versus t , for the LPF scheme. 
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The results for the IDF detection scheme show very similar characteristics to the 
results that are well known for PAM systems. This is to be expected, since a PAM 
detector usually employs an IDF, and it is the energy in each pulse that is measured. The 
LPF detection scheme does not measure the energy in each pulse, and therefore this 
method differs greatly from the IDF method. It is seen from the results that a LPF 


detector should be used to fully realise the advantages of PWM. 


CHAPTER V 
SUMMARY AND CONCLUSIONS 


In this thesis, a mathematical model of the noise present in multilevel PAM and 
PWM systems has been used to predict the bit error rate in a 4-level PAM system. In 
order to verify the accuracy of the theory, a 4.023 Mb/s, pseudo-random, 4-level PAM 
sequence was transmitted over a .5 km. multimode fiber and measurements were made 
of the error probability of the system as a function of the avalanche gain and optical 
power incident on the photodiode. 

The mathematical model makes use of the Gram-Charlier series to describe the 
PDF of the output voltage. The Gram-Charlier series takes into account the randomness 
of the avalanche gain, the amplifier's thermal noise and can include as well any generalised 
processing filter in the system. In addition to considering shot noise and thermal noise, 
the effect of intersymbol interference was also introduced into the Gram-—Charlier series 
model. 

Multilevel transmission is known to be advantageous compared to binary systems, 
since more information can be transmitted for a given channel bandwidth. The use of 
more than 4 levels, however, introduces a large amount of intersymbol interference. As a 
consequence, the analysis was restricted to 4-level PAM and PWM. 

Numerical examples have been provided to assess the performance of 4-level 
PAM and PWM systems. The system parameters that were used in these examples 
(Tables 3.1 and 4.1) represent realistic values. 

The optimal pulse amplitudes and thresholds were calculated for a range of error 
probabilities for a multilevel PAM system. The results obtained were given in Table 3.2. 
The pulse amplitudes and thresholds are seen, in this table, to be unequally spaced. To 
obtain the minimum error probability, the thresholds have to be chosen to lie closer to 
the lower level signals, rather than half-way between the signal levels. This results from 
the fact that the higher level optical signals have more shot noise than lower level signals. 

Graphs of error probability as a function of the average avalanche gain, <G>, of 
the APD and of the average optical power, Favg , reaching the receiver, are given in 
Figs. 3.4 and 3.5, respectively. It is seen in Fig. 3.4, that, for <G> less than 100, the error 


probability decreases as <G> increases; while for <G> larger than 100 the error 
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probability increases. This behaviour can be explained as follows: At low values of <G>, 
the thermal noise term dominates the overall noise. As <G> is increased, the output signal 
amplitude’ increases, but the thermal noise still dominates the shot noise. As a result of 
the thermal noise remaining constant while the signal increases, there is an overall 
increase in the SNR and, correspondingly, a decrease in the error probability. For <G> 
larger than 100, the shot noise starts dominating the overall noise. As <G> increases 
beyond 100, the shot-noise increases (the overall noise also increases) at a faster rate 
than the signal, thus resulting in an increase in the error probability. 

Fig. 3.5 shows that the error probability for a 4-level PAM system is highly 
dependent on the optical signal power reaching the detector. It is seen that a reduction of 
the optical power by only 3 db., causes the error probability to increase by four to five 
orders of magnitude. 

To verify the above theoretical results, an experimental optical-fiber 4-level PAM 
system was designed, constructed and tested at a transmission rate of 4.023 Mb/s. 
Measurements were made of the error probability as a function of <G> and Payg . and 
the results were plotted in Figs. 3.13 and 3.14. The theoretical calculations, based on the 
measured system parameters, were also plotted in Figs. 3.13 and 3.14. The theoretical 
and experimental results were seen to match fairly well. 

The experimental 4-level PAM system was seen to be thermal noise limited, as a 
consequence of which the probability density function of the output voltage was nearly 
Gaussian in shape. Only the first few terms of the Gram-Charlier series are neccessary 
to determine the bit error rates of this system. Thus the Gram-Charlier series model and 
the analytical model using the Gaussian approximation give almost identical results. To 
really utilise the advantages of the Gram—Charlier series one should consider a system 
with much less thermal noise. 

The analysis of the multilevel PWVM system proceeds in the same manner as for 
multileve! PAM. In PWM sytems, one is interested in preserving the timing content of the 
transmitted pulses. The pulse width is normally determined by threshold detecting the 
output pulse and measuring the width of the resulting regenerated pulse. Errors in 
determining the pulse width occur due to the noise present at the threshold levels. The 


timing error for the PWM system was expressed as a function of the noise amplitude. 
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Since the PDF of the noise amplitude can be characterised by means of a Gram-Charlier 
series, it is therefore also possible to characterise the timing error by means of a 
Gram-Charlier series. 

Two different PWM schemes, one employing a LPF and the other an IDF, were 
analysed. The optimum pulse-widths and threshold levels required to achieve a range of 
error—probabilities, for both the LPF and IDF schemes, have been tabulated (Table 4.2). 
Graphs of error probability as a function of <G> and P, were also plotted in Figs. 4.7 
and 4.8. From these figures one can observe the similarity in the characteristics of the 
PWM (IDF) system and the PAM system. This is to be expected, since the two schemes 
are basically the same. For both the PWM (IDF) scheme and the PAM scheme, the noise 
present during the whole pulse width is integrated and, depending on the value of the 
integrated noise, errors may or may not occur. The only difference between the PAM and 
PWM (IDF) scheme is that, whereas in the former, pulses of fixed width but different 
amplitudes are transmitted, in the latter, pulses of fixed amplitude but different widths 
are transmitted. The receiver section for both these systems is similar. 

The PWM (LPF) system, on the other hand, represents a true pulse width 
modulation system. The receiver section for the LPF scheme is entirely different from 
that used for the PAM system, since the pulse width is measured directly and is not 
determined indirectly from the amplitude of the signal output (as was done for the PWM 
(IDF) scheme). 

An intuitive understanding of the LPF scheme results from the analysis of equation 
(4.18). The value of K can be interpreted as the SNR of the LPF (PWM) system. It is seen 
that K is directly proportional to the signal amplitude, Vmjg, , but is inversely proportional 
to the standard deviation of the noise, a4 and the rise time of the received pulse, t,. 
Thus, for a LPF, as t, is decreased, the SNR increases and correspondingly the error 
probability decreases (see Fig. 4.9). One can, therefore, use very high bandwidth fibers to 
achieve low error rates even for very noisy systems. This is not possible for the IDF 
scheme because there is no such dependence of error probability on the rise time, t,. 

In conclusion, the LPF scheme can achieve very low error probabilities, provided 
that a wideband fiber is used [39]. The IDF scheme is, however, advantageous in systems 


where such wideband fibers are not available. 
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It is felt that the above results and the corresponding plots prove the usefulness 


of employing multilevel PAM and PWM transmission for fiber-optic communications. 


Though these systems have been analysed fairly comprehensively, there is scope for 


further research. The following improvements in the experiment are suggested: 


1. 


The pre-amplifier used in the experiment should be redesigned to provide much 
lower thermal noise. The Gram-Charlier series method could then be tested and 
compared to the other methods of calculating error rates. 

The data error analyser available for this experiment (HP 1645A) could operate only 
up to a frequency of 5 Mb/s. It is suggested that the BER measurements be 
repeated at higher bit rates, using a faster data error analyser. The advantages of 
multilevel transmission over regular binary transmission should become more 


noticeable at higher bit rates. 


Further research work could be undertaken in the following areas: 


i 


The analysis of the multilevel schemes was restricted to Gaussian shaped pulses. It 
would be interesting to note how the results would change if other pulse shapes, 
such as the raised-cosine were used. Impulse response measurements could be 
taken to determine the actual pulse-shape at the fiber output. 

The analysis of the PAM and PWM systems has not taken into account the effect of 
jitter, ie, the randomness of the timing signal. Jitter in the timing information, could 
strongly affect the error probability for commercial fiber-optic links. An analysis 
of jitter for PAM and PWM systems has been provided in [40]. The results obtained 
could be modified for multilevel systems. 

A comparison between multilevel PAM and PWM could be an interesting area for 
further research. An experimental PWM system could be designed and the results 
compared with those of the experimental PAM system, using the same system 


parameters. 


The author hopes that this thesis will provide a basis for further research work in the 


area of multilevel transmission for optical fiber Communications. 
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APPENDI®X | 
DERIVATION OF THE MOMENTS OF THE AVALANCHE GAIN <G >, FOR AN APD 

Personick [41], Mcintyre [42] and Dogliotti et a/ [26] have tackled the problem of 
deriving the moment generating function fon the avalanche gain of an APD. An overall 
view of the various approaches is presented here. Special emphasis is given to finding 
moments of the gain in a form suitable for the Gram-—Charlier series. 

The avalanche detection process essentially obeys Poisson statistics. Optically 
generated electron-hole pairs generate new electron-hole pairs through collision 
ionisation. These new electron-hole pairs can, in turn, generate additional pairs by the 
same mechanism. The avalanche process is depicted in Fig. A1.1. Charge carriers are 
created in the high-field region principally due to light pulses or thermal effects. Carriers 
generated outside the high field region can also drift into this region. In the high-field 
region, let the probability of holes suffering ionisation collision be B(x) and that for 
electrons be a(x), per unit length We define p, (nx) as the probability that n pairs 
ultimately result from the initially injected pair. The moment-generating function of the 


number of pairs, MIs), is therefore 
Mo(s,x) = x P;(Cn; x) ye" (A 1.1) 


The semi-invariant moment generating function (SIMGF), also called the cumulant, is 


defined as the natural logarithm of the moment-generating function and is given by: 


¥,(s,x) = In[M(sx)] = In[E, P,(n, x) e"] (A1.2) 


To derive the moment generating function, M (s,x), it is necessary to divide the high-field 
region into K intervals of width dX = W/K. These are labelled in the diagram (Fig. A1.1) as 
1,2,3,...,j..K. Assume that the initial electron-hole pair is injected into interval j. The 
probability density of a total number of pairs, n, resulting from the avalanche process 
(including the initial pair) is p, in.) with x taken as the centre of the interval j. 

The hole created in interval j moves to the left, towards x=0; the electron moves 
to the right, towards x=W. Both these holes and electrons can create new pairs in the 
other intervals in the high field region. The interval width, dX, is taken sufficiently small, 


such that either one or no new pairs are generated by the initial pair within this interval. If 
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Figure 1.1 Schematic of the avalanche region in an APD 
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the initial electron-hole pair generates a new pair in some other interval k, then this new 
pair will ultimately generate Ny, pairs, including itself due to the avalanche process. The 
total number of pairs generated is the sum of all the {Nx} produced in the different 
intervals. In order to get the SIMGF of the total number of pairs, we make use of the fact 
that the collision ionisation events are all independent (and therefore the Ny, are 
independent) and that the SIMGF of a sum of independent random variables is the sum of 
the individual SIMGF's. We denote the SIMGF of the total number of pairs by ie (s,x), 


i.e. 


¥,(s,x) =S + Ei w,(s) (A1.3) 


where: 
S = SIMGF of the initial pair and 
Pn, (s) = SIMGF of N. 
The SIMGF of N _ is obtained by noting that 


Yx(s) =Inf[E P,[N=zJe™ (A1.4) 
The probability that N, = Z is given by 
P,(Ny=z) = ydXP,(zx) for z>0 (A 1.5) 


The probability that Ny = O is givenby 1- dX 


where: 


a(x) if k>j 


of i 
B(x) if Kd 
and X = centre of interval k 


Fg (Z,x) = the probability that the total number of electron-hole pairs ultimately 


resulting in the avalanche process is Z. 


Equation (A 1.4) is simplified to 


os eee 


to we ht sv 26lgaigy imino raga | ta 
ee oF Gina accaenn saw ale aad end 
: . 
era) (aye “T+ tainty 
ey 
ae. ag Inithew arty Ne MMS Bee 
nis aaine 0 pa 
yartt ohineete ataini 
(hia ity (so KINA pnt = yt 5 
yal ravi ai S = gM wet eteiadong | 
ae ta O<s Tot (ng ythbe = fete | ; 


XB - GeouR eo 4g Wind 
eal sa | 
pad wi 


‘4 grns Yo tna © ki 
qa et 


= 


Verret Sieg sor=neyuels to wadrul ‘Rot ots 
ian 


¥,(s) = In[E, P.[N=zJ]e™ 
= In[ (1-ydX) oP P,(z, x) e"*ydx] 


In[ 1~ydX+ydXM,(s, x) ] 
In( 1) [ydX(e*e*-* -1) ] 


Taking the limit as dX becomes infinitely small and using (A1.3) and (A 1.6), 


¥.(s,x) = S+/o78(x' ) Letel* * 1) ] dx’ 
+ ff," a(x' )[ ete’ *)-1] dx’ 


Using Leibnitz's rule for the differentiation of integrals, one obtains 
0¥,(s, x) /Ox = [ B(x) a(x) ][e%""»-1] 
Equation (A1.8) can be solved to get 
¥,(s,x) = Inf{[1/{1—+Cexp{ fo™ [8(x") a(x") ]dx"}}]—1$ dx’ 
where 
es [eve(®-)-1] /etel*:%) 


Substituting (A 1.9) into (A1.7) the value of Vg (s,0) is obtained: 
¥,(s, 0) =$S + fp” a(zx')- 
{[1/{1-Cexp{ fo™ [8(x") a(x") ]dx"}}]—1} dx’ 


We make the assumption that at each point in the high-field region 


B(x) = ka(x) 


173) 


(A 1.6) 


(A 1.7) 


(A 1.8) 


(A 1.9) 


(A 1.10) 


(A 1.11) 


where k is known as the ionisation coefficient; ie, the ratio of the hole ionisation 


probability to the electron ionisation probability. 


One can then solve (A1.10) to obtain. 
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¥,(s, 0) = (s-d) + [1/(k-1)]- 
In{ e*)°/§ M.(s, 0) eT Ms, 0) -1] 3] 
= s4{1/(1-k) ]In] M,(s, 0) -e)*[ Ms, 0) -1]] (A1.12) 


where 
6 = fo" a(x)dx (the integrated ionisation 


coefficient for electrons) 
and M(s,0) = e%e*,° 
Differentiating (A1.12) results in 
@M,(s,0) /8s = M(s, 0) (k-1) /k]- 
tsi kK) TMs) 0) eaect cates 


(A 1.13) 


where MJ0,0) = 1 
Equation (A1.13) can be solved numerically. We can rewrite (A1.13) into a more 


suitable form by making use of new variables. 


M (s) = (k-1)M(s) /k[ 14 1/k) Me te6e) 9)] 
= pM(s) /{1-6[ M(s)e*] %} 


(A 1.14) 
where = (k-1) /k 
= 1 — (p/<G>) 
tS Te ea 
MO )iee al 
Differentiating (A1.14) we obtain: 
«3 Cinye 
12) M, 1 (A 1.15) 
where 
ag 2 
ag = 1-e-(&/p) 
(A 1.16) 
by = Be 


Using (A1.15) one can derive the following relation: 
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Differentiating (A1.15) iteratively and applying (A1.17) gives the I'th derivative of Mis) as 
22-1] 
yee les (1) \k. 
; bey *kg 


(A 1.18) 


where Ay, are functions of a,,a, and a;. 
2 


One can derive a relation between the coefficient Auay \- A k and Ay \-1 as: 


STs eee lel Ay a + Lia Cts Gls 2 ki Agee ae + (2-2)kA, o_o 
ee ae Saaeh oie Meas 
fe (A 1.19) 
A, : } a 0 
Finally, the moments of G can be written as: 
re 22-] 
Begs d°M(s Q 
a) SS 2 A <g> (A 1.20) 
iis" ie be K sk : 


where the relations, M(O} = 1 and M0) = <g>, have been used. 
Note that, for s=0, (a 15) gives an expression for <g’?> which agrees with 


U , 
Mcintyres excess noise factor [42]. 


F = [ <g*>/<g>7] = k<p> +.(2@-{ 17k) ) (i-k) (A1.21) 


The moments of the gain <g™, given in (A1.20), can be applied to (2.86) in order 


to obtain the cumulants of the received output voltage of a fiber-optic system. 
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APPENDIX II 
DERIVATION OF THE CHARACTERISTIC FUNCTION OF THE VOLTAGE AT THE 


OUTPUT OF AN OPTICAL-FIBER RECEIVER 
The aim of this appendix is to derive the characteristic function (CF) of the 
receiver output voltage. This is accomplished by first finding the CF of an electron-hole 
pair at the output of the detector. The effect of the avalanche gain on the CF is 


considered next. Finally the CF at the output of the filter is derived. 


FOF THE ELECTRON-HOLE PAIR AT TH TRUT OE THE DETECT 
Let g(t) be the shape of the optical pulse arriving at the input to the detector and 
let its width bet. The probability of n, electron-hole pairs being created in the 


sub-interval (t,, tp,, ), of length 4t,is given by the Poisson distribution: 


Px( Dy) = nek exp( —y,) /n,! (A2.1) 


where: py = ng(t,) AT/(hv) 
is the average number of charge carrier pairs in interval (ty. Ce ee): 
The electron-hole pairs created by the optical pulse, git), in different 
sub-intervals are independent of each other. The resulting output of the detector can be 
represented as a count vector, n = [n,,n,....., n,J, whose joint PDF, P, ("| is Obtained as a 


product of the individual PDFs, ie., 
P.(0,) = fl, Px( Dx) , (A2.2) 
The corresponding joint characteristic function, O,{jv), is given by: 


iy = = LJ D 
6,(jv) = fi 6,(jv) = tL, 2 0 exp(jvDy) Px( Mx) 
= fi, E exp(jvny) ek exp(—y) / ne 
= D; 
= ft (142% He u2/2) +--+) exp( He) 
(A2.3) 
But explx) = 1 + (x/ 1!) + (x2/2!) + (x3/3))+.. 


Therefore, (A2.3) reduces to: 
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@.(iv) = fl, exple] exp( 4) 
| = ft, exple"n — a] 
= exp{(e” -1) 2, iy, t (A2.4) 


F. AT THE FILTER OUTP 
The electron-hole pairs generated in the high field region of the p-n detector 
undergo avalanche gain multiplication. The CF of the final output electron-hole pairs 
after avalanche gain, in response to a single electron-hole pair, is given by (A1.12). This 


equation can be rewritten in the following form: 
In 6,(jv) = jv + (1-k)7* In[6@,(jv) (1-8) + a] (A25) 
where a = exp [ (k — 1) ] 


6 = fo" a(x)dx (wis the high field region widt h) 


where: 
o (x) is the collision ionisation probability for an electron. 
Assuming only electron ionisation events [ i.e. the factor k given in (A1.11) is neglected], 


the above equation can be rewritten as: 

8,(jv) = §1-G[ 1-exp(jv) ]}> (A2.6) 
where G = exp (& ) is the mean multiplication. 
For G> 10 , (A2.6) can be approximated as: 

6,(jv) = (1 — jvG) -1 (A2.7) 


Let P_(x) be the PDF of the number of charge carrier pairs after the random multiplication, 
in response to a single electron-hole pair. The PDF of the output in response to all the 


primaries in interval AT, is given by: 


P,(x) = Eo P,(0,) P(x) = (A2.8) 
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where: 
Bn 
p, x) = p (x)*p (x)*p (x)...%p (x) 
Let the characteristic function corresponding to p(x) be © (jv). Then 
O,(jv) = = £ Pal Xm) exp(jvXm) 
a ee a0 P,(n,) P,( x) 9 exp(jvx,,) 
Eo Pe(m) 2 Pi(xm)% exp(jvxn) 
Py(n,) exp[ n¥,(jv) ] 


Ey Pu(m) Oy (jv) = 5, 
Oy [ ¥.( jv) J 


(A2.9) 


where: 


YP liv) = In ©tivi 


For the time interval [0,t], the count vector, r, representing the random multiplication 


output, will have the joint characteristic function © (jv): 


é(jv) =f, al %(jv)] (A210) 
Since 
A o(iv) = exp [fexp(iv) 1} £, ud 


(A2.10) can be reduced to: 


exp [fexp(%(jv))-1} 2, ad 


6, (jv) 


But 


In@,( jv) 


%( jv) 
or 6,( jv) = exp ¥, (jv) 
6,(jv) = exp [{O(jv)) 1} By sa 
= exp [{0,(jv))-1} 49] 
(A2. 1 1) 


The characteristic function of the intersymbol interference and dark current electrons is 
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the same as (A2.11). Therefore, the APD output characteristic function [if we consider 


(p+q) interference Components plus dark current],is just a multiplication of terms like 


(AZ ly te, 


a(jv) = exp [{(Haldi-o) tat’ Ciudad}: 
i e.(iv) = 13) 


where u,! d,_,= average number of primary electron-hole pairs due 
to the bit under detection. 
U, = average number of primary electron-hole pairs due 
to the dark current 


Snlus d;) = the conditional intersymbol interference. 
4=-P 


(A2. 12) 


Equation (A2.13) is used in Chapter |! to derive the cumulant of the voltage at the output 


of the receiver. The characteristic furction of the Gaussian distributed thermal noise ts 


derived separately in Chapter Il. Since all the noise sources are independent, the overall 


characteristic function is the multiplication of the characteristic functions of the shot 


noise, intersymbo! interference and the thermal noise. 
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